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Abstract This paper considers the Fredholm determinant det(7 + XT of a Hankel integral 
operator on L 2 (0, oo) with kernel 4>(s+t+2x), where 4> is a matrix scattering function associated 
with a linear system (— A, B, C). The original contribution of the paper is a related operator 
R x such that det(J — R x ) = det(I — T x ) and —dR x /dx = AR X + R X A and an associated 
differential ring S of operators on the state space. The paper introduces two main classes of 
linear systems (—A, B, C) for Schrodinger's equation —if}" + utp = Xifj, namely 
Q\ . (i) (2, 2)-admissible linear linear systems which give scattering class potentials, with scat- 

tering function (j)(x) = Ce~ xA B; 
' (ii) periodic linear systems, which give periodic potentials as in Hill's equation. 

The paper analyses S for linear systems as in (i) and (ii), and the tau function is r(x) = 
det(I + R x ). 

(i) Here a Gelfand-Levitan equation relates </> and u(x) = —2^- logr(x), which is solved 
with linear systems as in inverse scattering. Any system of rational matrix differential equations 
gives rise to an integrable operator K as in Tracy and Widom's theory of matrix models. Under 
general conditions on existence of solutions, it is shown that there exist Hankel operators T$ 
and T^ with matrix symbols such that det(I + \xK) = det(J + //F$r^). The paper derives 
differential equations for r in terms of the singular points of the differential equation. This 

fSj 1 paper also introduces an admissible linear system with tau function which gives a solution of 

^ ■ Painleve's equation Pu. 

(ii) Consider Hill's equation with elliptic potential u. Then u is expressed as a quotient of 
tau functions from periodic linear systems. The general solution is a quotient of tau functions 
from periodic linear systems for all but finitely many complex eigenvalues A if and only if u 
is finite gap, hence has a hyperelliptic spectral curve. If the scattering data for the Gelfand- 
Levitan equation is a 2 x 2 symmetric matrix which is zero on the diagonal, then the solution is 
a symmetric matrix such that the on diagonal and off diagonal entries are related by differential 
equations that give the Miura transformation. 

The isospectral flows of Schrodinger's equation are given by potentials u(t, x) that evolve 
according to the Korteweg de Vries equation ut + u xxx — 6uu x = 0. Every hyperelliptic curve 
£ gives a solution for KdV which corresponds to rectilinear motion in the Jacobi variety of £. 
Extending Poppe's results, the paper develops a functional calculus for linear systems, thus 
producing solutions of the KdV equations. If T x has finite rank, or if A is invertible and e~ xA 
is a uniformly continuous periodic group, then the solutions are explicitly given in terms of 
matrices. 
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Key words: integrable systems, KdV, inverse scattering 



> 



i 



1. Introduction 



The motivation for this paper is from the theory of random matrices, and the scattering theory 
of differential equations with rational matrix coefficients. In Tracy and Widom's theory of 
matrix models [58], the basic data are a 2 x 2 rational differential equation and a curve. One 
starts with a system of differential equations 
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(1.1) 



with a, f3 and 7 rational functions, then one introduces a kernel 

f{x)g(y) - f(y)g(x) 



K(x,y) = 



x-y 



(1.2) 



which due to its special shape is known as an integrable operator. The other essential in- 
gredient of the theory is a prescribed union of intervals 7 = U^L [a2j-i, Ctej], so that K 
defines a trace class operator on L 2 (j); hence the Fredholm determinant det(7 — K) is defined, 
and one considers this as a function of the parameters dj. In particular, one can consider 
K : L 2 (0, 00) — > L 2 (0,oo) that is trace class and such that < K < I, so there exists a 
determinantal random point field on (0, 00), and det(J — KIt s oo \) is the probability that all 
random points are in (0, s). In applications to random matrix theory, the random points are 
eigenvalues of Hermitian matrices with random entries. 

Given an L 2 (0,oo) function </>, the Hankel integral operator F^ with symbol (f> can be 
defined on a suitable domain in L 2 (0, 00) by 



/>oo 

IV(*)= / <ttx + y)f(y)dy. 
Jo 



(1.3) 



When r^, belongs to the ideal c 1 of trace class operators on L 2 (0,oo), one can form the 
determinants det(/ + nT^) and the eigenvalues of G c 1 satisfy multiplicity conditions which 
are stated in [44, 48]. More generally, one can introduce <t>( x ){y) = 4>{ x + 2y) and consider 



t(x;h) = det(/ + /J> (a0 ) 



(1.4) 



as a function of x > and fi G C. In this paper, we analyse t(x, /i) by the methods of linear 
systems. Tau functions of this form arise in scattering theory of Schrodinger's equation. There 
are also applications to integrable operators, since for significant cases of (1.2), such as the Airy 
kernel or Bessel kernel [58, 59], there exists a Hankel integral operator such that F 2 = K; 
hence one can describe det(J — K) in terms of r(x,yu). In [9] we showed how one can realise 
r<£ by means of linear systems. In the present paper, we take linear systems as the starting 
point and show how general properties of the linear system are reflected in the r functions and 
systems of differential equations so produced. 

Definition (Linear system) Let H be a complex Hilbert space, known as the state space, 
and B(H) the space of bounded linear operators on H. Let (e~ tA ) t >o be a Co semigroup of 
bounded linear operators on H such that ||e~ M || < M for alH > and some M < 00. Let 
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T>(A) be the domain of the generator —A so that V{A) is itself a Hilbert space for the graph 
norm H^Hf^) = ||£Hh + ||-<4.£||hj an d let A^ be the adjoint of A. Let H be a complex separable 
Hilbert space which serves as the input and output spaces; let B : Hq — > H and C : H — >■ Hq 
be bounded linear operators. The linear system (—A, B, C) is 

^ = -AX + BU 
at 

Y = CX, X(0) = 0; (1.5) 

so (f>(x) = Ce~ xA B is a bounded operator function on Hq, and the corresponding Hankel 
operator is F^ on L 2 ((0, oo); H ), where T 4> f(x) = / °° <j>{x + y)f(y) dy. 

Definition (Admissible linear system). Let (—A,B,C) be a linear system as above; suppose 
furthermore that the observability operator ©o : L 2 ((0, oo); Hq) — >■ H is bounded, where 

6 /= / e- sAt C^f(s)ds; (1.6) 

suppose that the controllability operator E : L 2 ((0, oo); Hq) — >■ H is also bounded, where 

/>oo 

H /= / e- sA Bf(s)ds. (1.7) 

(i) Then (—A, 5, C) is an admissible linear system. 

(ii) Suppose furthermore that ©o and So belong to the ideal c 2 of Hilbert-Schmidt oper- 
ators. Then we say that (—A, B, C) is (2, 2)-admissible. 

In [9, Proposition 2.4] we showed that for any (2, 2) admissible linear system, the operator 



/•oo 

R x 



poo 

= / e- tA BCe~ tA dt (1.8) 

J X 

is trace class, and the Fredholm determinant satisfies 

det(/ + Ai2 x ) = det(/ + Ar^ (a . ) ) (x>0,AeC). (1.9) 

Whereas R x does not have a direct interpretation in control theory, the notation suggests that 
R x has many of the properties of a resolvent operator, as we justify in Lemma 2.2 below. In 
examples of interest in scattering theory one can calculate det(J + XR X ) more easily than the 
Hankel determinant directly [31, 48]. The operator R x has additional properties that originate 
from Lyapunov's equation, and which make it easier to deal with than r^ (x) . 

Definition (Lyapunov equation). Let —A be the generator of a Cq semigroup on H and let 
R : (0, oo) — >■ ~B(H) be a differentiable function. The Lyapunov equation is 



dz 

with initial condition 



dR 

" ! '' z =AR Z +R Z A (z > 0) (1.10) 



AR + RqA = BC. (1.11) 
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The definition differs slightly from the equations from [44, 48] that define the observability 
and controllability Gramians. In this paper we take (1.10) as the starting point, instead 
of r</, {:i;) ). In section 2 we solve (1.10) for some (2,2) admissible linear system, then use 
R x to construct solutions to the associated Gelfand-Levitan equation; thus R x is a tool in 
the scattering problem. The following definition of u is motivated by scattering theory for 
Schrodinger's equation —ip" + uip = Aifi in L 2 (R). See [21] 

Definition (Potential). For each (2, 2) admissible system with Hq = C, the potential is 

U{X) = lo S det ( 7 + r ^))- ( L12 ) 

Theorem 1.1 (i) Suppose that (—A, B, C) is a (2, 2) admissible linear system with A bounded. 
Then there exists a solution R x to (1.10) and (1-11) such that t(x) = det(7 + R x ) is entire. 

(ii) Alternatively, suppose that (—A, B, C) is a linear system with input and output space 
H, and (e~ xA ) is a uniformly continuous and it-periodic group on H. Suppose that there 
exists a trace class operator E on H such that AE + EA = BC. Then there exists a solution 
to (1.10) and (1.11) such that t{x) = det(J + Rx) is entire and it-periodic. 

(Hi) In either case u is meromorphic on C. 

Part (i) is proved in section 2, while (ii) is proved in section 7. In [11] we introduced 
examples of periodic linear systems as in (ii), and here develop a systematic theory which 
shares some common elements of scattering theory from case (i). 

The fundamental idea of [41] is to realise Hankel operators with balanced linear systems; 
we refine this idea by working with admissible linear systems, so that we can define deter- 
minants and hence the r function. In section 2, we solve the Gelfand-Levitan equation by 
means of the operator R x and recover u from <f>, as in inverse scattering. Tau functions natu- 
rally arise in parametrized families. In particular, the tau functions tq(x) = det(J + R x ) and 
Ti(x) = det(7 — Rx) are linked by the Miura transformation, as we discuss in section 2. 

To realise integrable operators as in (1.2), we need to work with products of Hankel op- 
erators. Poppe [41, 49, 50] proved some remarkable product formulas involving products and 
traces of Hankel integral operators and applied them to scattering theory, and his work moti- 
vated some of the results of this paper. In section 3, we introduce a functional calculus which 
encompasses Poppe's ideas, but uses R x and operators on the state space H of (—A,B,C). 
We suppose that (e~ tA ) defines a holomorphic semigroup and we can introduce a domain £1 on 
which det(7 + R z ) is holomorphic and nowhere zero, so I + R z has a bounded inverse F z . We 
introduce a differential ring S of holomorphic functions Q — > B(H), which contains A, BC, R z 
and F z , so that we can solve (1.10) and (1.11) inside S. If we can choose S to be a right 
Noetherian ring, then we say that (—A, B, C) is finitely generated. Given S, we introduce a 
space of functions B and the linear map [ . J : S — >■ B such that 

[P\ =-^tmce(p(F x -I)). (1.13) 
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We identify a subring A of S such that the range of [ . J restricted to A is a differential ring 
|_AJ of functions which contains u(x). In these terms, the scattering transform is 



cf)( x ) = Ce- xA B < — > u{x) = -4[A\ . (1.14) 

Thus |_ • J linearizes the determinant formula in (1.12). 

Gelfand and Dikii [26] considered the ring Ao = C[u,u' , -u" , . . .] of complex polynomials 
in u and its derivatives. They showed that if u satisfies the stationary higher order KdV 
equations (3.15), then — /" + uf = \f is integrable by quadratures on a spectral curve, which 
is a hyperelliptic Riemann surface £ of finite genus. Such u are known as finite gap or algebro 
geometric potentials since + u has a spectrum in L 2 (R) that consists of intervals known 
as bands, separated by finitely many gaps. Then A is a Noetherian ring; see [15, 53]. The 
ring |_AJ is analogous to Ao in the particular examples that we analyse in subsequent sections. 

In section 4 we show that if A is a finite matrix with eigenvalues Xj such that 9fJAj > 0, 
then (— A, B, C) is finitely generated. We also recover some determinant formulas from the 
theory of solitons. 

In section 5 we show how to realise kernels of the form (1.2) from linear systems by means 
of products of Hankel operators with matricial symbols. The system of differential equations 
(1.1) depends upon the poles of a, j3 and 7, hence these are natural parameters for the solution 
space. We recall how Schlesinger's equations [51, 22] arises in this context, and compare various 
notions of tau functions by the partial differential equations that they satisfy. 

Krichever and Novikov [38] considered 

[JL- Uj ,L]=B J L (1.15) 

where Uj are matrix functions and Bj are differential operators, a relation which is similar 
to (1.13). They formulated the notion of an algebo-geometric system. In particular, this 
applies to finite gap Schrodinger equations, where the spectral parameter may be chosen to 
be a meromorphic function on a hyperelliptic Riemann surface. We recall that a compact 
Riemann surface £ is hyperelliptic if and only if there exists a meromorphic function on £ that 
has precisely two poles. In this case, there is a two-sheeted cover £ — >■ P 1 with 2g + 2 branch 
points, where g is the genus of £. The elliptic case has g = 1. 

Our main application of differential rings for integrable operators is in section 6, concerning 
the Airy kernel 

^ 2 Ai(x)Ai'(y)-Ai'(x)Auj/)) 

—y ' (1-16) 

where <j>{x) = Ai(x) is the Airy function. This is is a universal example in random matrix 
theory [59], and in particular i^Oc) = det(J — } /4) is the cumulative distribution function 
of the Tracy- Widom distribution associated with the soft spectral edge of the Gaussian unitary 
ensemble. We recover Ablowitz and Segur's result of [2] that v(x) = — 2(logi*2(x)) // satisfies the 
Painleve's second transcendental differential equation Pu and v(x) decays rapidly as x — >■ 00. 
Airy's equation is the first of a sequence of differential equations </>( 2 ^(x) = xcj)(x), and we show 
that the corresponding 77 functions satisfy the Painleve hierarchy of differential equations. 
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The T£ are linked by the string equations, which in this case is given by the Gelfand-Dikii 
polynomials [26, 27]. These facts are well known, but our method proof also gives a natural 
probabilistic meaning to the T£. 

In discussing Hill's equation, Ercolani and McKean [21] observe that the notions of Jacobi 
variety and theta functions can be extended to the case of infinitely many spectral gaps, 
whereas the notion of a multiplier curve is somewhat tenuous. Likewise we can introduce 
tau functions via determinants of linear systems in cases where there is no related algebraic 
curve of finite genus. The spectral class of a potential is invariant under flows associated with 
the Korteweg de Vries equation ut + u xxx — 6uu x = 0, which belongs to a hierarchy of partial 
differential equations which are themselves associated with flows u(0, x) u(t, x) on the space 
of potentials. In section 8, we consider a class of differential rings that are especially suited for 
describing solutions of KdV. When U (t) is a family of unitary operators on H, the tau function 
and potential of (—A,U(t)B,CU(t)) evolves with respect to t. In section 7, we introduce a 
unitary group U (£) so that </>, u and [ . J itself evolve with respect to t as in the KdV flow, and 
thus we linearize the the KdV flow on functions of rational character, and produce solutions of 
the higher order KdV equations. As is standard in this context, we consider a family of t(x; t) 
functions parametrized by t G T°°, the infinite-dimensional torus. 

If u is a finite gap potential for Hill's equation, then the spectral curve is hyperelliptic 
and has a finite-dimensional complex torus X as its Jacobi variety, thus the corresponding tau 
function can be expressed as the restriction of a theta function to a straight line in the tangent 
space of X by results of Its and Matveev [26, 27]. In section 8 we formulate a sufficient 
condition for the tau function of a linear system to be algebraic, in this sense, in terms of 
the Kadomtsev-Petviashvili equations [62] . Soliton solutions of KP occur for spectral curves 
that are rational curves in the plane that have only regular double points. The term elliptic 
solitons refers to functions of rational character on the torus, namely elliptic functions, which 
give solutions of the KdV equation. 

A significant advantage of R x and S is their application to periodic linear systems, which 
seem to lie outside the scope of [41, 49]. In section 9, we introduce linear systems (—A, B, C) 
such that A is an invertible operator that commutes with BC, and e xA is a uniformly contin- 
uous periodic group and the A, B, C are block diagonal matrices. Thus we introduce periodic 
linear systems with potentials that are either rational trigonometric functions on the complex 
cylinder C/nZ or elliptic functions on the complex torus C/tt(Z + iirZ) as in section 11, and 
show that these have analogous properties. 

The table below summarizes the functions that we produce from explicit linear systems 
in sections 5,6,9 and 11. Here g is the genus of the spectral curve, p is Weierstrass's elliptic 
function, 6\ is Jacobi's theta function [42], u in the fifth column satisfies Pji from [22]. 



equation 
Schrodinger 



u G LAJ 
scattering 



r G L 



£ 

R -f [0,oo) 



Painleve 
Hill 
Lame 
soliton 



g(g + l)cosech 2 x 



finite gap 
-9(9 + l)p 



Pi 1 



Tracy-Widom F 2 

e 

(sinhx)f(9 +1 )/ 2 



hyperelliptic 
{- 5 ,...,-l}U[0,oo) 
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In section 10, we introduce the family of linear systems T,\ = (—A, (XI+A)(XI—A)~ 1 B, C) 
for A in the resolvent set of A, and the corresponding tau function t\(x); then we introduce 
the Baker- Akhiezer function iIjba{x, A) = e Xx T\(x)/r(x); here x is the state variable and A a 
spectral parameter; see [6, 35,37] for related notions of ipBA- We say that (Sa)a is a Picard 
family of linear systems if ar i — >■ iPba(x, A) is meromorphic for all but finitely many A G C. This 
term is introduced by analogy with the terminology of Gesztesy and Weikard [28, Theorem 
1.1], who define a meromorphic potential u to be Picard if — /" + uf = Xf has a meromorphic 
general solution for all but finitely many A G C. Indeed, u is finite gap if it satisfies the 
stationary KdV equations as in [26, 27]. We therefore consider a family of linear systems 
T,\(t), with common A : H — >■ H, and constant input and output spaces, where t = (ti, ti, ■ ■ ■) 
is a sequence of real parameters and A is a spectral parameter. Then £,\(t) has a potential 
u\(x;t) with poles depending upon (A, t); thus the dynamics of the system is reflected in the 
pole divisor of the potentials, as we describe in section 9. For various linear systems, we 
introduce a compact Riemann surface £ and a meromorphic function A : £ — >■ P 1 such that 
A i — ^ iPba(x, A) is meromorphic, except possibly at finitely many points. 

Our most complete results are for elliptic potentials, as in section 11. We obtain a char- 
acterization of the elliptic potentials that are finite gap in terms linear systems. All elliptic 
potentials can be realised as quotients of r functions from periodic linear systems, however, the 
general solution of Hill's equation can be expressed as a quotient of r functions from periodic 
or Gaussian linear systems only if the potential is finite gap. This complements results of 
Gesztesy and Weikard from [28]. In [33], Kamvissis recovered the determinant formula (1.12) 
for scattering potentials as a limiting case of the Its-Matveev formula for periodic finite-gap 
potentials as the period tends to infinity. In this paper, we show that the periodic linear system 
which gives rise to Lame's equation also gives a kernel which degenerates to Carleman's kernel 
under thermodynamic and high density limits. 



2 Solving Lyapunov's equation and the Gelfand-Levitan equation 

We begin with simple existence result, showing how linear systems in continuous time give 
rise to Hankel matrices. Subsequent results will introduce stronger hypotheses to ensure the 
existence of Fredholm determinants. 

Proposition 2.1 Suppose that H is a separable Hilbert space, and that 

(i) C : H — >■ C and B : C — >■ H are bounded linear operators; 

(ii) A is a densely defined linear operator in H; 

(in) A is accretive, so &{Af, f) > for all f G T>(A); 

(iv) XI + A is invertible for some A > 0. 
Then (e~ )t>o is a Co contraction semigroup on H, so (x) (s) = Ce-( 2x+s ^ A B is bounded 
and continuous on (0, oo); the cogenerator V = (A — I) (A + J)~ 1 satisfies \\V\\ < 1 as an 
operator on H, and there is a unitary equivalence between r^ (x) on L 2 (0, oo) and the Hankel 
matrix on £ 2 (N U {0}) that is given by 

CO 



V2Ce- 2xA V n+m (I + A)~ 1 B 
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Proof. By the Lumer-Phillips theorem [20], —A generates a Co contraction semigroup. Di- 
rectly from the definition (iii) of an accretive operator and hypothesis (iv), one proves that 

11*11 < !■ 

We introduce the Laguerre polynomials of order zero L { n\s) = (n\)- 1 e s (d/ds) n s n e- s and 
then the functions h n (s) = \/2e~ s Ln\2s), so that (h n )™ =0 gives a complete orthonormal basis 
of L 2 (0, oo). By integrating by parts, one can verify that 



Jo y/2n\ Jo ds n 



s n e- s ds 



ds n 

= V2Ce- 2xA (A-I) n (A + I)- n - 1 B. (2.2) 

Peller [45, p. 233] shows that F^ is unitarily equivalent to the Hankel matrix under the unitary 
correspondence (h n )^ =0 {e.j)JL , where (e^) is the standard orthonormal basis of £ 2 . 

□ 

We introduce Lyapunov's equation, and the existence of solutions for suitable (—A, B, C). The 
solution R x is defined by a formula suggested by Heinz's theorem [8, Theorem 9.2] and has 
properties analogous to the resolvent operator of a semigroup. 

For < 6 < 7r, we introduce the sector fie = {z G C \ {0} : | argz\ < 9}. 

Lemma 2.2 Let (—A,B,C) be a linear system such that He - * " 4 )! < 1 for some to > 0, and 
that B and C are Hilbert-Schmidt operators on Hq such that ||.B||ifs||C||ifs < 1- 

(i) Then (—A, B, C) is (2, 2) -admissible, so the trace class operators 

[>(X> 

R x = e~ tA BCe- tA dt (x > 0) (2.3) 

J X 

give the solution to (1-10) for x > that satisfies (1.11), and the solution to (1-11) is unique. 

(ii) Suppose further that there exist M and tt/2 < 6 < ir such that XI + A is invertible 
and satisfies \X\ \\ (XI + A)^ 1 \\ < M for all X G fV Then R z extends to a holomorphic function 
which satisfies (1.11) on Qo_ n / 2 , and R z — >■ as z — >■ oo in for all < e < 6 — tt/2. 

Proof, (i) Since BC G c 1 , the integrand of (2.3) takes values in c 1 and is weakly continuous, 
hence strongly measurable, by Pettis's theorem. By considering the spectral radius, Engel and 
Nagel [20] show that there exist 6 > and M s > such that \\e- tA \\ < M 5 e- St for all t > 0; 
hence (2.3) converges as a Bochner-Lebesgue integral with 



/>oo 

i4|| c i < / Ml\\BC\\ c ie- 2St dt 



' X 

r2 



M 2 

< ^ll^lks||C||^e^. (2.4) 
Furthermore, A is a closed operator and satisfies 

,4 / e- tA BCe- tA dt+ C \- tA BCe~ tA dtA = C -^e- tA BCe~ tA dt 



dt 

= e- xA BCe~ xA - e- TA BCe~ TA 

-). e- xA BCe~ xA (2.5) 
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as T ->• oo; so AR X + i? x ^4 = e^^BCe -2 ^ for all x > 0. We deduce that x •->■ R x is a 
differentiable function from (0, oo) to c 1 and that the modified Lyapunov equation (1.8) holds. 

Now suppose that AR + R A = BC and AW + W A = BC, and consider V = R - W . 
Then for £,77 G H, we have 

j t (V Q e- tA t, e- tA \) H = ((V A + AV Q )e~ tA ^ e~ tA \) H = 0; (2.6) 

hence (Voe~ tA ^, e^ tA \)u is constant, and by the hypothesis on A, we have 

(Vbe~ tj4 £, e~ tA] rj)H — >■ as t — >■ 00. Hence (Vo£, 77) fj = 0, and so Vq = 0, and i?o is unique. 

See [41, p. 261] for a similar argument. 

(ii) By classical results of Hille, e~ zA defines an analytic semigroup of 0, e _ n / 2 , so we can 
define R z = e~ zA Rue~ zA and obtain an analytic solution to Lyapunov's equation. For all 
< e < 8 - 7r/2, there exists M' £ such that ||e~ zA || < M' £ for all z G Sle- £ —K/2, so R z is 
bounded on this sector and by (2.4), R z — >■ as z —> 00 in 

□ 

We now introduce and interpret the operator R x and the companion r function in the 
context of Proposition 2.2. The operators R x satisfy a linear differential equation and the r 
function is given by an expectation of a random variable with respect to a reference measure 
on state space. Let H be a separable Hilbert space, and let R x be a trace class operator on H. 
Suppose that E is a positive and trace class operator on H, and that L x is a Hilbert-Schmidt 
operator on H such that I + L x is invertible, L X E = EL X and (I + L X )(I + R X )E{I + R x ) = E. 

Proposition 2.3. (i) The operator E x = (I + R X )E(I + R x ) is also positive and trace class 
on H. 

(ii) The operators E x and E determine equivalent Gaussian probability measures j x and 
7 on H with mean zero, so that 7^ and 7 have the same null sets. 

(in) The Jacobian of the linear transformation between 7^ and 7 is t(x) = det(7 + R x ). 

Proof, (i) Let I - K x = (I + L^ 1 , so that K x is a Hilbert-Schmidt operator on H, with 
one not in the spectrum of K x , and K X E = EK\. Evidently E x = (I — K X )E is positive and 
trace class. 

(ii) Hence there exist Gaussian probability measures 7 and 7^ on H with mean zero and 
with characteristic functions 

/ e i ^ w) lx {dw) = e- {E ^ )/2 (£eH) (2.7) 

J H 

and likewise with 7 for j x and E for E x . To check that the measures have the same class of 
null sets, we introduce the Hilbert space / D(E^ 1 ^ 2 ) = {£ G H : £ = E 1 / 2 ^, for some ry G H} 
with the norm ||£||^ B _ 1/2 j = (E~ l £,t;). One can introduce a complete orthonormal basis 

for H consisting of eigenvectors of E, and hence produce an orthogonal basis for "D^E^ 1 ^ 2 ); 
then one checks that K and EK^ E^ 1 give Hilbert-Schmidt operators on T>(E^ 1 ^ 2 ). By the 
general theory of Gaussian measures, the cr-algebra of null set of 7 is determined by the pair 
V{E~ 1 / 2 C H, and likewise the a-algebra of null sets of 7^ is determined by V(E X 1 ^ 2 ) C H, 
and we have shown that these pairs of Hilbert spaces are equivalent. See [39]. 
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(iii) We can write (I+R x )£ = r\ and then (E x £, = (En, 77), so this linear transformation 
takes j x to 7, and the Jacobian in det(7 + R x ). 

□ 

Definition (i) A bounded linear operator Y on Hilbert space is said to be block Hankel if there 
exists 1 < m < 00 such that Y it is unitarily equivalent to the block matrix [Aj+k~2]j°k=i 011 
£ 2 (C m ) where Aj is a m x m complex matrix. 

(ii) Let (—A, B, C) be a (2, 2) admissible linear system with input and output space Hq, 
where the dimension of Hq over C is m < 00. Then m is the number of outputs of the system, 
and systems with finite m > 1 are known as MIMO for multiple input, multiple output, and 
give rise to block Hankel operators with = Ce~ xA B; see [48]. 

(iii) The Gelfand-Levitan integral equation is 

/>oo 

T(x,y) + ${x + y) + fi / T(x, z)${z + y) dz = (0 < x < y) (2.8) 

J X 

where T(x, y) and $(x + y) are m x m matrices with scalar entries. 
Proposition 2.4 (i) In the notation of Lemma 2.2, there exists x > such that 

T M (x, y) = -Ce- xA (I + nR^e^B (2.9) 

satisfies the integral equation (2.8) for x < x < y and < 1). 

(ii) The determinant satisfies det(I + iiR x ) = det(J + nY^, (x) ) and 

yutraceT M (x,x) = ^- logdet(J + fiR x ). (2.10) 

Proof, (i) We choose x so large that e 5x ° > Ms/25, then by (2.4), we have I^HI-R^H < 1 for 
x > xq, so I + [iR x is invertible. Substituting into the integral equation, we obtain 

Ce -(x+y)A B _ Ce -xAy + ^R x )-^ e -vA B 

/>oo 

- nCe~ xA (I + /iik) -1 / e-^SCe"^ dze~ yA B 

J X 

= Ce- {x+y)A B - Ce- xA (I + nR^e'^B - ^iCe~ xA (I + [iR^R^^B 
= 0. (2.11) 

(ii) As in (1.6), the operator : L 2 (0,oo) — >■ H is Hilbert-Schmidt; likewise S x : 
L 2 (0, 00) ->■ i? is Hilbert-Schmidt; so (-A,B,C) is (2, 2)-admissible. Hence r$ = Oj,^ 
and i? x = S X ©J. are trace class and 

det(J + /ii2x) = det(J + ^6*) = det(J + //O^) = det(J + fiT^). (2.12) 
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Correcting a typographic error in [9, p. 324], we rearrange terms and calculate the derivative 



fiTpfax) = -^trace (Ce _xyl (J + fxR^e-^B^j 
= -^trace(7 + fxR x )~ 1 e~ xA BCe~ xA 

(It t-> \ — 1 dR x 

(1 + fiR x * 



dx 



= — trace log(7 + uR x )- 
dx 



(2.13) 



This identity is proved for < 1 and extends by analytic continuation to the maximal domain 

of Tfj,(x, x). 

□ 

Now we show how to calculate the determinants of Hankel products in terms of the 
Gelfand-Levitan equation. Changing to a more symmetrical notation, we suppose that 
(—A\,B\,C\) and (—A 2 ,B 2 ,C 2 ) are (2,2) admissible systems with state spaces Hi and H2 
and output space that realise (f>\ and (f> 2 . First, let Rjk ■ H k — > Hj for j, k = 1, 2 be the 
operators 

/>oo 

R jk (x)= e~ tA * BjC k e- tAk dt {x > 0). (2.14) 

J X 



For the first result, we introduce that state space H = 
H = C 2xN and A : H ^ H , B : H ^ H and C : 77 — » i7 by 

A = 

so that 



#1 
H 2 



~A 1 


" 


, B = 


'B 1 


" 


, c = 


" 


c 2 





A 2 _ 





B 2j 








and the output space 



(2.15) 



<E>(x) = 



fa(x) 
Mx) 



(2.16) 



Proposition 2.5 (i) There exists 5 > such that for all fx E C such that \fx\ < 5, I 
fx 2 R 2 i(x)Ri 2 (x) has an inverse G x and 



T(x,y) = 



fxC 2 e- xA *G x R 2X (x)e-y A ^B x -C 2 e- xM G x e~ vM B 2 
-C ie - xA i (7 + fx 2 R 12 {x)G x R 21 (x))e-y A i B 1 fxC ie - xA ^ R 12 {x)G x e~y A * B 2 

(2.17) 

satisfies (2.8) for all x > xq from some xo > 0. 
(ii) The determinants satisfy 



det(7 - fx 2 R 12 {x)R 21 {x)) = det(7 + Mr* ((B) ). 

and their logarithmic derivatives are 

d 



dx 



logdet(7-/x F^ :(x) r ( ^ 1 = ^traceT(x,x). 
11 



(2.18) 



(2.19) 



(Hi) In particular, with A 2 = A\, B 2 = eC\ and C 2 = B\ and e = ±1, the identities hold 
with 4> 2 {x) = £0i (x)^ ; hence if e = 1, then T$ is self-adjoint; whereas if e = — 1, then F$ is 
skew. 

Proof, (i) It is easy to check that = Ce~ xA B. Likewise, we can compute 



Rx — 



- tA BCe~ tA dt 



R 12 (x) 
R 2 i(x) 



(2.20) 



which is a trace class operator on H since both (—Ai,Bi,C\) and (— A 2 ,B 2 ,C 2 ) are (2,2)- 
admissible. For x so large that |//| 2 ||i?i2(x)|| ||i?2i(a;)|| < 1, we can form the operator G x = 
(7 — /u 2 i?2i-Ri2) _1 and hence compute 



F x = 



I nRi2{x) 
ldR 21 (x) I 



I + u. 2 R 12 (x)G x R 21 (x) -u,R 12 (x)G a 



-fxG x R 21 (x) 



G 3 



(2.21) 



Then we compute T(x,y) = —Ce~ xA F x e~ yA B and obtain the matrix from (2.9). One then 
checks, as in Lemma 2.2, that T satisfies the integral equation (2.8). 

(ii) We introduce the observability operators Q x : 7 2 ((0, oo); C 2xAr ) — >■ 77 by 



Ox 


7" 




' 


e 2 " 




7" 




_9 _ 




_0i 







_9_ 




and the controllability operators 


^x 


7 2 ((0 


oo 


);C 2 




7" 




^2 


" 




7" 




^x 
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Si 









Ce-*^C 2 g(t)dt' 

r e -tAi c t mdt 



Ce- tA *B 2 f(t)dt 
Ce-^B l9 (t)dt 



(2.22) 



such that 



as operators on 77, and 



o ~ 2 e\ 












R21 




Rl2 






e}-! 







°2,(x) 







.ejs 2 

as operators on 7 2 ((0,oo); C 2xN ). Now from the determinant identity 

det(7 + fi,GlE x ) = det(7 + /iS^e* ) 

we deduce 

det(J - ^T^r^) = det(7 - u 2 Ri 2 (x) R 2 i(x)) . 
The function i2 x is differentiable with respect to x, so by Lemma 2.1, we can compute 

— logdet(7-/i 2 r^ 2(x) r^ 1(:c) ) = — logdet(I + fiR*) 

= /itraceT(x, x). 



(2.23) 
(2.24) 

(2.25) 

(2.26) 
(2.27) 



(2.28) 
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(iii) We have fafa) = C 1 e~ xAl B 1 and fa(x) = C 2 e- xM B 2 = eB\e- xA ^C\. 

□ 

Remarks (i) In applications to physics and probability, it is important to identify those 
such that r<£ = and r(x) = det(J + r^ (a;) ) is real for real x. 

Dyson considered scattering theory with even potentials u G (R, R) with no bound 
states, which can be realised from linear systems with m = 1 and self-adjoint F^, as in [9, p. 
324]. By increasing rank to m = 2, we can allow finitely many bound states. 

The case of m = 2 includes models associated with the Zakharov-Shabat system [62]. 
Within this group, there are scattering functions of the form 

<f> + iip~ 

<\> — tip o 

such that (j) an d ip are real, so det(7 — T) = det(7 — r^ +i ^r^ +i ^,); see section 6 of [9]. 
Definition Suppose that t\ and To are meromorphic functions that satisfy 

v = (d/dx) log(ri/r ), w = (d/dx) log(r ri), u = -2(d 2 /dx 2 ) logr 



$ = 



and that 



r^n - 2t^t[ + t t[' = 0. 



(2.29) 



Then u is the Miura transform of v; see [30]. 

Furthermore, ro and r\ satisfying (2.29) are actually holomorphic. Moreover, if tq has 
a zero of order a at zo and t\ has a zero of order f3 at zo, then local analysis shows that 
(a - f3) 2 = a + p. 

In the following result, we show how products and quotients of r functions can be linked 
by the Gelfand-Levitan equation for 2x2 matrices. If the scattering data for the Gelfand- 
Levitan equation is a 2 x 2 symmetric matrix which is zero on the diagonal, then the solution 
is a symmetric matrix such that the entries on the diagonal and the entries off the diagonal 
are related by differential equations that give the Miura transformation. 

Theorem 2.6 Let {—A, B, C) be a (2, 2) -admissible linear system with input and output 
spaces C, and let 4>{x) = Ce~ xA B. 

(i) Then there exists S > such that for all \i G C such that |//| < 5, the integral equation 
() with 



T(x,y) = 



has a solution such that 



W(x,y) V(x,y) 
V(x,y) W(x,y) 





(x + y) 



4>{x + y) 




(2.30) 



W(x,x) = logdet(7 
dx 2/i 



(2.31) 



v(x ,)-A± log det ( / + ^)) 

V[X ' X) dx2^ l0g det(I-^T^ ) )- 



(2.32) 
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and 

±-j-W{x,x) = -V(x,x) 2 . (2.33) 

(ii) If 4> is real, then there exist complex numbers {±Zj(x) : j = 1,2,...}, mutually 
distinct, such that 

det(J J Wgy i ) _ " / l + ^-frU 
det(/-^ w ) Lj\l-u.z j (x))- 

(iii) If fi = 1/2, then u{x) = -2(d 2 /dx 2 ) logdet(J - 2" 1 r^ (;c) ) is the Mi ura transform of 
V(x, x). 

Proof, (i) This is a special case of Proposition 2.5 in which 4>i = fa = 4>i so we can take 
A\ = A 2 = A, Bi = B 2 = B and C\ = C 2 = C, so we write R\ 2 = R 2 i = R- Then we obtain 
the required solution with diagonal terms 

W(x, y) = ixCe- xA {I - /i 2 R 2 ) -1 R x e~ yA B (2.35) 

and off-diagonal terms 

V(x, y) = -Ce- xA (I - /j 2 R 2 )~ l e~ yA B. (2.36) 
Hence the off-diagonal terms satisfy 

V(x, x) = -trace((/ - ix 2 R 2 x y 1 e~ xA BCe~ xA ) 

= itrace(((7 - pR^ 1 + (/ + yRx)' 1 ) ^Rx) 

= ^^(logdet(/ + u.R x )- log det(I-u.R x )y (2.37) 

Next we observe that 

V(x, x) 2 = Ce- xA (I - /j 2 R 2 x )-\AR x + R X A)(I - /i 2 Rly 1 e~ xA B , 

and by repeatedly using Lyapunov's equation (1.10), one can check find 4-W(x, x) from (2.35); 
thus one obtains (2.33) Finally, we can express the Fredholm determinants in terms of the 
Hankel operator r^ (x) . 

(ii) We suppress the dependence of the terms upon x. Since F^ is trace class, the 
spectrum consists of together with non-zero eigenvalues Xj with multiplicity ^({Aj}), and 
for a self-adjoint the algebraic and geometric multiplicity are equal. Then the function 
fi i— >■ det(J + /jTfj,)/ det(J — /iF^) is meromorphic on C, and the formal difference of its zeros 
and poles on {/i 6 C : |//| < p) may be represented as a divisor ■|<p( z/ ({Aj})^-i/A ;/ ~~ 
u ({^j})$i/\j)- Evidently the coefficients satisfy Y2j-\i/\ j \< p (. l/ (.{^j}) ~~ = 0> where 
^({Aj}) — ^({—Xj}) belongs to { — 1, 0, 1} by a theorem of Megretskii, Peller and Treil [44]. 

(iii) One can check that (2.29) holds if and only if u = v' + v 2 and v 2 = —w'. 

□ 
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The differential equation 




/ 



9 



,W = 



v 



; k 
k —v 



(2.38) 



is equivalent to — /" + (v' + v 2 )f = k 2 f and kg/ f = f'/f — v, where u = v' + v 2 is the Miura 
transform of v. 

The main application is in the context of the (modified) Korteweg-de Vries equation. The 
mKdV equation is associated with a pair of r functions, namely tq and t\ which satisfy (2.29) 
and give rise to a solution v. Then the Miura transform u satisfies KdV, and is associated with 
To only, so the relationship between KdV and mKdV is not symmetrical. When describing 
matrix models for quantum field theory [30], one says that the partition function for KdV 
is the square of one r functions; whereas the partition function for mKdV is the product of 
two r functions which are linked by (2.37) as in []. In section 6, we obtain scaling solutions 
of mKdV from scaling solutions of the second Painleve equation by constructing a matricial 
linear system as in Theorem 2.6. 

In the next section we show to interpret such calculations systematically. 

3 The state ring associated with an admissible linear system 

In the established theory of linear systems, one uses the Laplace transform to move from the 
state variable, and then one considers factorization in rings of transfer functions as in [23] . In 
this paper, we prefer to work with differential rings of operators on the state space so as to 
integrate various differential equations related to Schrodinger's equation. We introduce these 
state rings in this section, and develop a calculus for R x which is the counterpart of Poppe's 
functional calculus for Hankel operators. As we see in subsequent sections, our theory of state 
rings has wider scope for generalization. 

Definition (Differential rings) Let H and K be separable complex Hilbert spaces, let B(iJ) 
be the ring of bounded linear operators on H . For xq, X\ G R let S be a subring of 
C°°((0, oo); B(H)); that is we suppose that each T G S is a differentiable function of x G (0, oo) 
as we indicate by writing T x ; we suppose further that dT x /dx G S, and that (d/dx)(ST) = 
{dS/dx)T + S(dT/dx). Then S is a differential ring with the subring {S G S : dS/dx = 0} of 
constants. When / G S, we identify 91 with 9 to simplify notation. 

Definition (State ring of a linear system) Let (—A, B, C) be a linear system such that A G 
B(iJ). Suppose that: 

(i) S is a differential subring of C°°((0, oo); B(H)); 

(ii) /, A and BC are constant elements of S; 

(iii) e~ xA , R x and F x = (I + R x )~ l belong to S. 
Then S is a state ring for (—A, B, C) on [xq, X\). 

(iv) Moreover, if S is left Noetherian as a ring, then we say that (—A, B, C) is finitely 
generated. 

Remarks. When A is algebraic, we can use simple functional calculus to help construct the 
differential ring. We use this technique in sections 4 and 9. 
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Lemma 3.1 Suppose that (—A, B, C) is a linear system with bounded A and that R x gives 
a solution of Lyapunov's equation (1.10) such that I + R x is invertihle for x > with inverse 
F x . Then the free associative algebra S generated by I, Rq, A, Fq, e~ xA , R x and F x is a state 
ring for (—A, £?, C) on (0, oo). 

Proof. First we note that BC = AR + R A belongs to S, as required. We also note that 
(d/dx)e~ xA = —Ae~ xA and that Lyapunov's equation (1.10) gives 

^(1 + R^- 1 = (I + R x y\AR x + R X A)(I + R x y\ (3.1) 
dx 

which implies 

dF 

—^ = AF X + F X A-2F X AF X . (3.2) 
dx 

with the initial condition 

AF + F A - 2F AF = F BCF . (3.3) 

Hence S is a differential ring. 

□ 

Definition (Complex differential rings and state rings) Let 17 be a domain in C and 'M.q(X) 
the meromorphic functions from 17 to some complex Banach algebra X. If S as above is also 
a subring of M.q(X), then we use the standard complex derivative d/dx and say that S is a 
complex state ring for (—A, B, C) on Q. 

In the context of Lemma 2.2(h), one can introduce xq > such that ||-Ri || < 1, then use 
the sector ft = {x = xq + z E C : $tz{Af, f) > 0,V/ G T>(A)} so that F x is holomorphic on 
17. In section 7, we work with periodic meromorphic functions and replace 17 by the complex 
cylinder C/ttZ. In section 11, we work with doubly periodic and meromorphic functions, so 
we replace 17 by T = C/A, where A is a lattice.) 

Definition (Brackets) Given a state ring S for (—A, B, C), and let B be any differential ring 
of functions from (0, oo) — > B(K). Let 

A = span c {A n \ A ni F x A n2 . . . F x A nr : nj G N}. (3.4) 

Then let [X, Y] = XY - YX be the usual Lie bracket on S, and let [J : S — >• B be the linear 
map 

[FJ = Ce- xA F x YF x e~ xA B (Y G S). (3.5) 



Lemma 3.2 (i) Then A defines a differential subring of S. 

(ii) The range |_SJ is a differential ring with derivative d/dx, and has [A\ as a differential 
subring. 

(Hi) Suppose that S is contained in c 1 . Then 

trace[Pj = trace(P(AF x + F X A - 2F X AF X )). 
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In particular, suppose that the input and output spaces are C. Then 



[P\ = trace(X(AF x + F X A - 2F X AF X )). 

Proof, (i) The bracket operation satisfies 

[P\ [Q\ = [P{AF X + F X A - 2F X AF X )Q\ , (3.6) 



dP 



— [P\ = [A(I - 2F X )P + — + P(I- 2F X )A\ . (3.7) 

We can multiply elements in S by concatenating words and taking linear combinations. Since 
all words in A begin and end with A, we obtain words of the required form, hence A is a 
subring. To differentiate a word in A we add words in which we successively replace each F x 
by AF X + F x A — 2F X AF X , giving a linear combination of words of the required form. 

(ii) As in (i), the operations are well defined in the sense that [P\ [Q\ and (d/dx) |_PJ are 
images of elements of A for all P, Q G A. Evidently the proposed multiplication is associative 
and distributive over addition. Using (3.6) and (3.7), one checks that Leibniz's rule holds in 
the form 

d (WL«J) = (|:W)L«J + W^L«J). (3-8) 



dx 



(iii) We have 



trace(P— -) =tracePF x e- xA BCe- xA F x 

= tvace(Ce- xA F x PF x e- xA B), (3.9) 

and F' x = AF X + F X A - 2F X AF X , whence the result. 

Finally, observe that when C has range in the scalars, we can remove the trace and write 



trace 



dP \ 

P^p: j = Ce- xA F x PF x e- xA B = [P\ . (3.10) 



□ 

Definition (Liouvillian extension) Let K be a field of complex functions with differential d, 
and adjoin a complex function h to K where either: 

(LI) h = J g for some g G K, so that dh = g; 

(L2) h = exp J g for some g EK; 

(L3) h is algebraic over K. 
Then K(/i) is a Liouvillian extension of K as in [60]. More generally, a field L is a Liouvillian 
extension of K if there exist differential fields Fj such that K = Fo C Fi C . . . C F n = L, 
and each Fj arises from Fj_i by applying (LI), (L2), or (L3). 

Theorem 3.3 Let (—A, B, C) he a linear system as in Lemma 2.2, and suppose furthermore 
that A is hounded and H = C. 

(i) Then (—A, B, C) has a complex state ring S on C on which R z is unique. 

(ii) The map [.J : S ->■ M C (C) satisfies <f>(2x) = [P^ 2 J and u(x) = -4[A\. 
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(iii) The ranges |_SJ and [A\ are differential rings. The held of fractions K of |_AJ is a 
differential held, and t(x) = 1/ detF x is entire and belongs to a Liouvillian extension L of K. 

(iv) C(u, v! , . . . , u( k ~^) is a differential subheld of K, if and only if = r(u, . . . , u^ fc_1 ^) 
for some rational function r. 

Proof, (i) Mainly this follows from Lemma 2.1 and Proposition 2.4. By Riesz's theory of 
compact operators, the F x = (I + R x )~ defines a meromorphic operator valued function on 
C. Hence we can select S to be the subring of meromorphic functions from fl to B(H) generated 
by /, A, BC, R x , e~ xA and F x . On {x : R x + R x > —21}, the function F x is holomorphic and 
satisfies F' x = FA + AF — 2FAF. 

(ii) Evidently 

[F~ 2 J = Ce- 2xA B = cf)(2x), (3.11) 

while we can write (1.12) as u(x) = —2(d/dx) 2 \ogdet(I + R x ) where (d/dx) logdet(J + R x ) = 
[F~ 1 \ and differentiate using (3.2). 

(iii) From the definition of R x , we have AR X + R X A = e~ xA BCe~ xA , and hence 

F x e- xA BCe~ xA F x = AF X + F X A - 2F X AF X , (3.12) 

which implies 

[P\ [Q\ = Ce- xA F x PF x e- xA BCe- xA F x QF x e- xA B 

= Ce- xA F x P(AF x + F X A - 2F x AF x )QF x e^ xA B 

= [P(AF X + F X A - 2F X AF X )Q\ . (3.13) 

Moreover, the first and last terms in \_P\ have derivatives 

^-Ce- xA F x = Ce- xA F x A(I - 2F X ), -^—F x e~ xA B = (I - 2F x )AF x e~ xA B, (3.14) 
dx dx 

which implies (3.7). Hence by Lemma 3.2(h), the image of |_-J is a differential ring. 

Now |_AJ is a subring of Mc(C) and hence is an integral domain with a field of fractions 
K. We have 2(d/dx) 2 logdet F x = u(x) G K, so we can recover det F x by integration and 
exponential integration. By (2.3) and Morera's theorem, R x is an entire c^-valued function, 
hence det (7 + R x ) is entire. 

(iv) By (iii), u and all its derivatives belong to K. Evidently C(u, . . . ,u'' k ^ 1 ^) is a dif- 
ferential field if and only if it is closed under differentiation, or equivalently, equals some 
rational expression in u, . . . ,u ( > k ~ 1 \ 

□ 

Definition (Stationary KdV hierarchy) (i) Let gi = — (l/4)w. Then the KdV recursion formula 
is 

d d d d 3 

A dx~ 9m+l ^ = 8gi ^dx~ 9m ^ + 8 dx~( 9l ^ 9m ( X ^ + dx^ 9m ^' ^ 3 ' 15 ^ 
The solutions may depend upon constants of integration; if the constants of integration are 
chosen all to be zero, so that 52 = (3/16)n 2 — (l/16)u" etc, then the g m give the homogeneous 
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KdV hierarchy, and g m is a differential polynomial in u, known as a Gelfand-Dikii polynomial. 
In this case, the differential equations g m = are known as Novikov's equations; see [26, 27]. 

(ii) If u satisfies g m = for all m greater than or equal to some mo, then u satisfies the 
KdV hierarchy and is said to be an algebro-geometric (finite gap) potential. 

Suppose that u is finite gap. Then the KdV hierarchy has a solution in the Noetherian 
differential ring C[u, u',u", . . .] of differential polynomials in u. More generally, we will ob- 
tain commutative Noetherian state rings in which we can solve the KdV hierarchy, using the 
following result. 

Proposition 3.4 Let S be a commutative differential ring given by a linear system (—A, B, C), 
and suppose that [P\ = Ce- xA F x PF x e- xA B. Then 

4 |^ 2m+3 J = ^~ 3 [A\ +8(^[A\)[A^\ +W[A^\(^[A\). (3.16) 

Proof. We can obtain the following identities by repeatedly using the basic calculus rules 

— \A\ = \2A 2m+4 - 2A 2m+3 FA - 2AFA 2m+3 \ ; (3.17) 
ox 

\A\ A \A 2m+1 \ = \ 2AFA 2m+3 + 2A 2 FA 2m+2 - 2AFA 2m+2 FA 
ox 

- AAF AF A 2m+2 - 2A 2 FAFA 2m+l - 2AFA 2 FA 2m+1 

- 2A 2 FA 2m+1 FA + 4AFAFAFA 2m+1 + 4AFAFA 2m+1 FA\ ; (3.18) 

( — [A^j [A 2m+1 \ = [2A 3 FA 2m+1 + 2A 2 FA 2m+2 - 4A 2 FAFA 2m+1 - 4AFA 2 FA 2rn+1 

- AAFAFA 2m+2 + 8AFAFAFA 2m+1 \ ; (3.19) 



| A 2m+i | = 1 8A 2m+4 _ 2 4A 2m+3 FA - 24AF <A 2m+3 - 24A 2m+2 FA 2 - 24A 2 FA 2m+2 

dx 3 

+ A8A 2m+2 FAFA + 48AFA 2m+2 FA + 48AFAFA 2m+2 - 8A 2m+l FA 3 

- 8A 3 FA 2m+1 + 24A 2m+1 FA 2 FA + 24A 2m+1 FAFA 2 + 24A 2 FA 2m+1 FA 
+ 2AAF A 2m+l F A 2 + 2AA 2 FAFA 2m+l + 24AF A 2 F A 2rn+1 

- A8A 2m+l FAFAFA - 48AFA 2m+1 FAFA - 48AFAFA 2m+1 FA 

- 48 AF AF AF A 2m+l \ (3.20) 

We now use the assumption that S is commutative to simplify terms, and obtain the stated 
result by cancelling. 

□ 

Remarks 3.5 (i) To deal with calculations as in Theorem 2.6, one can adjoin G = (I — R) -1 , 
which satisfies G' = GA + AG — 2GAG. (ii) Airault, McKean and Moser [3] consider the cases 
of Theorem 3.3(iv) given by v!" = \2uv! for u rational, trigonometric and elliptic. 
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(ii) Poppe [49, 50] introduced a linear functional [.] on Fredholm kernels K(x,y) on 
L 2 (0,oo) by \K~\ = K(0, 0). In particular, let K,G,H,L be integral operators on L 2 (0, oo) 
that have smooth kernels of compact support, let V = T,p (x) have kernel 0(s + t + 2x), let 
T' = j^Y and G = T^, (x) be another Hankel operator; then the trace satisfies 

fTl = - -^-tracer (3.21) 

ax 

\FKG] = ---^-traceTKG (3.22) 
2 dx 

\(I + T)- l T] = -trace((/ + r)- 1 r'), (3.23) 

\KT] \GL\ = -\\K{T'G + TG')L] , (3.24) 

where (3.24) is known as the product formula. The easiest way to prove these is to observe 
that T'G + YG' is the integral operator with kernel — 20(^ (s)^^ (t), which has rank one, as 
in (6.?) below. These ideas were subsequently revived by McKean [41]. 

(iii) Mulase [45] considers differential rings over C that are also closed under (LI) and 
(L2); an important example is the Noetherian ring C[[x]] of formal complex power series. 
However, C[[x]] does not contain functions with poles. Krichever [35] considered an algebraic 
curve with a preferred point po, and functions that are holomorphic except for poles at po- 
Note that {f(z) = ^2 < kL^ n dkZ k ;n G N; a& G C} is a Noetherian differential ring, but it is 
not closed under (LI) or (L2). So we prefer to start in a smaller ring and then control the 
extensions that are formed by making quadratures. 

(iv) Let (S, d) be a differential ring with / G S, and let G = {I + £7=-oo : aj G S} 
be the formal Volterra group with the integration operator We wish to consider the 
operator T as in the Gelfand-Levitan equation and select S so that I + T E G, where I + T 
is the analogue of the wave operator of scattering theory. To obtain an S with a relatively 
simple form, we suppress the input and output spaces of (—A, S, C) and deal with operators 
on H. Then we use a special bracket operation |_ . .J to return to scalar- valued functions and 
we identify a differential field that includes r. 

In the next three sections, we give significant examples of differential rings associated with 
linear systems. 

4 Finite matrix models 

In this section, we are concerned with complex differential rings for linear systems (—A, B, C) 
that have finite dimensional state spaces. While we seek to realise S by functional calculus, 
we do not assume commutativity of A and BC, and we do not assume that e~ xA is stable. 

Hypotheses. Throughout this section, we let A be a n x n complex matrix with eigenvalues \j 
with geometric multiplicity rij such that Xj+Xk / for all j and k; if all the eigenvalues are geo- 
metrically simple, then let K = C(e~ Al *, . . . , e~ Xnt ); otherwise, let K = C(e~ Al *, . . . , e~ Xnt ,t). 
Also, let B = (bj) G C nxl and C = ( Cj ) G C lxn . 

The following result extends a special case of the Sylvester-Rosenblum theorem [8] . 
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Lemma 4.1 Let S = C{I, A, BC}. Then there exists R G S such that RqA + ARq = -BC, 
and the equations (1.10) and (1.11) have a unique solution. 

Proof. Let E be a chain of circles that go once round each Xj in the positive sense and have 
all the points —Xk in their exterior. Then by [8], the matrix 

R = J- [ (A + XI)- 1 BC(A- XI)- 1 dX (4.1) 
2m Jz 

gives the unique solution to the equation (1.11). Furthermore, by the Cayley-Hamilton theo- 
rem, (A =F XI)- 1 is a polynomial in A, A, I and det(A =p XI)- 1 for all A on 7; hence i?o belongs 
to the algebra So- 

The function R x = e- xA R^e- xA is entire and of exponential growth, and gives a solution 
of (1.11) and (1.10). Since R x is of exponential growth, it has a Laplace transform R(s) = 
Jo°° e ~ sx Rxdx which satisfies sR(s) + AR(s) + R(s)A = R , and for all s > 2\\A\\ the solution 
is unique and may be expressed as 

/ioo 7 \ 

((A + 8/2)1 + Ay'Rod-X + 8/2)1 + A)' 1 — . (4.2) 
-ioo ^ m 

Hence R x is the unique solution of (1.10) and (1.11). 

□ 

Theorem 4.2 Let R x = e- xA R e- xA ; then let S = K{I, A, BC}. 

(i) Then (—A,B,C) is finitely generated since S is a left Noetherian ring with respect to 
the standard multiplications. 

(ii) The linear map [ . J : S -)> H c : [P\ = Ce- xA F x PF x e- xA B satisfies 0(2z) = [F- 2 \ 
and u(x) = — 4 . Also, r, r/r A G K. 

Proof (i) The complex algebra generated by /, BC and A is finite-dimensional and hence left 
Noetherian; so by Hilbert's basis theorem, S as a subalgebra of M n (K) is also Noetherian; see 
[15, p. 106]. Observe that (XI - A) (XI + A)- 1 G S for all -A in the resolvent set of A. 

By the Riesz functional calculus, we can introduce a sum of cycles going round each Aj 
once in the positive sense, so that 




hence there exist complex polynomials pj and qj, and integers mj > such that 

n 

e- tA = ^ qj (t)e-^ Pj (A), (4-4) 
i=i 

where qj(t) is constant if the corresponding eigenvalue is simple. Hence R x G S, and likewise 
all the entries of R x belong to S. Moreover, for any B G C nxl and C G C lxn , there exist 
constants aj and polynomials qj such that 

n 

<j>(x) = Ce- xA B = <\!<lM)< X - r - (4.5) 
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Now introduce the minors Uj G K of i" + R x such that 

det(/xJ - (/ + R x )) = u n + an^ix)^- 1 + ... + ai (xV + {-l) n e(x), (4.6) 

and recall that by the Cayley-Hamilton theorem 

(/ + R x ) ((J + it,)"" 1 + (7 n _i(a;)(/ + R x ) n - 2 + ... + <n(x)ij + (-l) n £(x)I = (4.7) 

so F x belongs to S. Hence S is a complex differential ring for (—A,B,C). By the usual 
expansion of the determinant, r G K. 

(ii) This follows as in Theorem 3.3. Observe also that (f> and u belong to K, and all 
elements of K are meromorphic on C. 

□ 

Lemma 4.3 (The Cauchy determinant formula) Let x r and y s be complex numbers such that 
x r y s 7^ 1. Then 



det 



1 - Xjy k 



j,k=i 



Yil<j<k<n( X j X k)Y\.l<m<p<niym Up) 
ril<7\s<n(l ~~ X rUs) 



(4.8) 



Proposition 4.4 Suppose that B = (bj)™ =1 G C n , C = (cj)™ =1 G C lxn and A is the n x n 
diagonal matrix with simple eigenvalues Xj such that Xj + A*; / for all j = 1, . . . , n. 
(i) Then R x gives rise to the determinant 



bee 

det (J + iiR x ) =l + nJ2 iCjC 



-2A,-a: 



2A,- 



bibrr,CkC n e ( X j+ X k+\n + >~p)x 



(j,k),(m,p):j^m;k^p 



(Aj + A m )(Afc + A p ) 



j = l l< i<fc <n ^ + Ak) 



(4.9) 



Proof, (i) The proof is by induction on n. There is an expansion 



det 



Ojk + 



u.b j c k e-( x j +x ^ x 



Xj + Xk 



= V ^ CT det 

<rC{l,...,n} 



Aj + Afc J j.fciEcr 



(4.10) 



in which each subset a of {1, . . . ,n} of order jjcr, contributes a minor indexed by j, k E a. 
Letting x r = A r and y r = — 1/X r in the Cauchy determinant formula, we obtain the identity 



det 



bjC k e 



-XiX — XhX 



Xj + Afc J j,k£a 



2A 



JA A, - X k 

3 j,k<Ecr:j^k 3 ^ k 



n 



(4.11) 
□ 
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Remarks 4.5 (1) The preceding results of this section apply in particular when A is a finite 
matrix such that all the eigenvalues have > 0. 

(2) Kronecker's theorem asserts that a bounded Hankel integral operator has finite rank if 
and only if the transfer function is a rational function with all its poles in {z G C : ^Rz < 0}. 
Such rational functions are known as stable. In [23] , the authors consider factorization of the 
transfer function in M nxn (C(A)) and the subring of stable matrix rational functions. Their 
results describe the properties of S rather than S itself. 

As we mentioned in Lemma 2.1, Hankel integral operators are unitarily equivalent to 
infinite Hankel matrices, which naturally compress to sequences of finite Hankel matrices. We 
now state some results concerning sequences of r functions from finite Hankel matrices. 

Definition Let V G B(H) satisfy ||V|| < 1, and let C : H -)• C and B : C -)> H be bounded 
linear operators. Then = Ax n + Bu n and y n = Cx n define a linear system in the discrete 
time variable n G N, and there is an associated Hankel matrix V = [CV^ + 2 B]j° k=1 . The 
principal minors r n = det[CV J+ 2 -B]™ fc=1 give the corresponding tau sequence (1,ti,T2, . • .)• 

We now consider how some finite tau determinants can arise from linear systems on an 
infinite-dimensional Hilbert space. Suppose that {^j)JL Q is a complex sequence such that the 
Hankel matrix V = [7j+fc-2]£fc=i defines a bounded linear operator on H. Then by a theorem 
of Young [44], there exist bounded linear operators V : H — >■ H , B : C — >■ H and C : H — >■ C 
such that 7j = CV^B for j = 0, 1, . . .. Then we introduce To = 1, t\{x) = Ce~ xV B and 

r n (x) = det[Ce- xV V^ +k - 2 B]l k=1 . (4.12) 
Proposition 4.6 (i) Toda's equation holds 

= 3-^. (4 ,3) 

(ii) If (L;d/dx) is a differential Geld that contains t\(x), then r n (x) belongs to L for 
n= 1,2,.... 

Proof, (i) We can express the tau functions in terms of minors of derivatives 

r n ( X )=det\ \ xj+ t 2 X) ]., , (4.14) 



j,k=i 



Then by a formula of Darboux [24, 61, 32], the second derivative of logr n (x) is a multiple of 
the exponential of the second difference of the sequence \ogr n (x). 
(ii) This follows from the recurrence relation in (i). 

□ 

Suppose further that V = V' is bounded on H and B = C'; then there exists a scalar 
spectral measure p on the spectrum a of V, and a measurable family (iJ^)^ ecr of Hilbert spaces 
such that H = J® H^p(d^) and Vv = £v for all v G H^. Determinants of finite Hankel matrices 
appear in random matrix theory via the eigenvalue distributions of random Hermitian matrices, 
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especially in generalized unitary invariant ensembles. Let M^(C) = {X G M n (C) : X = X^}, 
and for X G M^(C), let Ai < . . . < A n be the eigenvalues, listed according to multiplicity. 

Now let p be an absolutely continuous probability measure on [—1,1], and for t = 
(ti,t2,...) G £°°(R), let pt(dx) = exp(J2J^ 1 tjX 2: >~ 1 )p(dx). Then the sequence of moments 
p m (x;t) = J^ l X m p t ( y dX) gives rise to a Hankel matrix T(x;t) = [fJ>j+k-2(x',t)]j° k=1 ; then the 
principal minors 

r n (x;t) = det\jjL j+k - 2 (x;t)]l k=1 (n = l,2,...) (4.15) 

give tau functions; let r = 1. Let [z\ = (2z~ 2j+1 /(2j - l))^i, so that t ^ t + [z\ gives 
the (odd terms in) the Miwa-Sato shift. In this context t\ is a time variable; z is a spectral 
parameter; n is the number of rows in the matrix, and n \— > n + 1 is referred to as a Backlund 
transformation as in [1, 31]. 

Proposition 4.7 There exist a bounded and positive measure v n j(dX) on M^(C) such that 

(i) v Uj t(dX) is invariant under the action of the unitary group on M^(C) by X (->■ UXU^ ; 

(ii) the ratio of tau functions satisfies 

T n (x;t+[z\) _ / { X6MA(C):A w <s} det ((^ + X )( zI ~ ^P^ Rt^X) 



T n{x]t) /{l£Mj(C):A„<i} v n,t{dX) 

Proof, (i) Let A : M^(C) — >■ R be the map that associates to X its eigenvalues Ai < 
A2 < • • • < A n , listed according to multiplicity. Note that the set {X G M^(C) : X n < x} 
is invariant with respect to unitary conjugation. Then there exists a bounded and positive 
measure u njt (dX) on M n (C) that is invariant under unitary conjugation and such that A 
induces 

^n >t (rfA)= Yl (Xj - X k f ptidXx) . . . p t (dX n ) (4.17) 

l<j<k<n 

on R n , where the product is essentially the Jacobian of the transformation. One can express 
this in terms of Vandermonde's determinant, and deduce that 

f a n , t (dX) = [ detfAj- 1 ] 2 p t (dX 1 )p t (dX 2 ) . . . Pt (dX n ) 

J(-l,x) n J(-l,x) n 



eX 

= det 



X^- z p t {dX) 

=T n (x;t). (4.18) 

(ii) We now make use of the special form of p t , and reduce the left-hand side to 

„ 00 

T n (x;t)= / exp(Vt J trace(X 2 ^ 1 ))i/ n>0 (rfX) (4.19) 

J{xeMZ((cy.\ n <x} j=1 

and likewise 

r n (x;t+[z\) (4.20) 
„ 00 00 

= / exp(2 V^ 1 " 2i trace(X 2 ^ 1 )/(2j - 1) + Vt j trace(X 2 ^ 1 ))i/ n0 (dX) 

J{XeMK{C):\ n <x} j = 1 j = 1 
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so we obtain the stated result by forming the ratio of these. 



□ 



In the same way, the determinants of finite Toeplitz determinants appear in random matrix 
theory from random unitary matrices, especially in the circular ensembles. 

5 Tracy— Widom kernels and Schlesinger's differential equations 

In [19], Deift, Its and Zhou observe that the correlation functions for many exactly solvable 
models admit determinant representations involving integrable operators, although only a few 
of them are of convolution operators. In random matrix theory, one often encounters integrable 
kernels that are the products of Hankel integral operators on I/ 2 (0, oo); see [58, 59] and section 

6 below for examples. The most significant examples in random matrix theory are associated 
with linear systems of differential equations with rational matrix coefficients, and on this 
section, we introduce (2,2) admissible linear systems with vectorial input and output spaces, so 
that we can realise integrable operators as products of Hankel operators with matrix symbols. 
Then we consider the properties of the tau functions that arise from these Hankel operators, 
and derive Schlesinger's system of differential equations. 

Definition (Integrable operators) [19] An integrable kernel has the form 



K(x,y) = 



J2l=ifj(x)gj{y) 



x 



(5.1) 



where fj,gj are continuous and bounded functions on (0, oo), and we suppose further that 
Y^j=i fj( x )9j( x ) = 0, so K is nonsingular on x = y. 

In particular, consider the system of differential equations 



dx 



= tt(x) 



n(x) 



7 a 
a (3 



J 



-1 

1 



(5.2) 



with a, /3 and 7 rational functions. Then, as in Tracy and Widom's theory of matrix models 
[58, 59], we introduce the kernel 



f{x + 2z)g(y + 2z) - f(y + 2z)g(x + 2z) 



x 



(5.3) 



and L (2) by (I - L (z) )(I + K [z) ) = I. 

Theorem 5.1 Suppose that a, f3 and 7 are proper rational functions with n poles of order less 
than or equal to p, for some p G N, and all poles are in C \ [0, 00); suppose that f,g<EL 2 (0, 00) 
are solutions of (5.2) and that f(x),g(x) —> as x — >■ 00. 

(i) Then there exist Hilbert-Schmidt Hankel operators F$ and with 2np 2 x 2np 2 matrix 
symbols $ and ^ such that 

det(I + XK (z) ) =det(/ + Ar* w r* w ). (5.4) 

(ii) There exists xq such that L( 2 ) is a bounded integrable operator for all z > xq. 
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(Hi) Suppose further that e 2ex f{x) — )■ and e 2£X g(x) — >• as x — >■ oo for some £ > 0. Then 
$ and \I/ are realised by (2, 2) admissible linear systems. 

Proof, (i) We can write 



n Pk 



(5.5) 



where the i?o and for £ = 1, . . . ,pk with Pk < p and = 1, . . . , n are symmetric 2x2 
matrices and the poles dj lie in C \ [0, oo). From the differential equation, we have 



V dx dy/ 



x-y 





~m~ 




~f(y)~ 


y 


_g(x) _ 




Ay) . 



(5.6) 



n Pk £ 



= -£££ 



k=l £=1 v=0 



E, 



kl 



(x - a k ) 



t-v 



g{x) 



(y - ak) 



u+1 



f(y) 
g(y) 



where we have used the real inner product. Noting that Ek,t has rank less than or equal 
to two, let N = 2np 2 and introduce scalar- valued functions <f)j(x) and ij)j(y) such that the 
previous sum equals — X^jLi ( t > j{ x )ipj(y), and since the poles are off (0, oo), we can ensure that 
J °° x(\4>j(x)\ 2 + \ijjj(x)\ 2 )dx is finite, so <j>j and ipj give the symbols of Hilbert-Schmidt Hankel 
operators on L 2 (0, oo). Then one verifies the identity 



f(x)g(y) - f(y)g( 



x 



x-y 



N 



^2<f>j(x + s)ipj(s + y)ds; 

3=1 



(5.7) 



indeed by the preceding calculation, the difference between the two sides of (5.7) is a function 
of x + y, which goes to zero as x — >■ oo or y — >■ oo in any way. and hence vanishes identically. 
Finally, we build the JVxJV matrices 



<3>(x) 



M x ) <f> 2 (x) ... 4>n{x) 
... 















My) o ... o 
My) o ... o 

4> N (y) o ... oJ 



(5.8) 



so that T$ and are Hilbert-Schmidt matrix operators, and with 4>j,{ z ){x) = (f>j(x + 2z) etc 
we have 



det 



N 

(I + \K (Z) ) =det(/ + Aj]r^. (z) r^ ) ) =det(/ + AT« w r» w ). 



(5.9) 



Hence 



H(x) = 



A$(x) 
-*(&) 



has det(I + XK (z) ) = det(/ + r s ). 
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(ii) We can define Li z \ = Kr z \(I + Kt z \) 1 for all z such that < 1. Now let 8 be 



any derivation on the bounded linear operators on I/ 2 (0, oo), and observe that 

SL=(I + K)-\5K){I + K)' 1 . 



(5.10) 



In particular, with Mh(x) = xh(x) for h G L 2 (0,oo), the derivation 5K = MK — KM 
is represented by the finite-rank kernel f(x)g(y) — f{y)g(x) which vanishes on the diagonal 
x = y; hence ML — LM is also a finite rank kernel which vanishes on the diagonal. In short, 
we obtain L from the kernel 



F(x)G(y) - F(y)G(x) 
x-y 



where 



F 




[OH 




G 






-K)- l g_ 



(5.11) 



Moreover, SK = [d/dx, K] is the finite rank integral operator that is represented by the kernel 
(5.6), so 5L is also finite rank. 

(iii) Given that / and g are of exponential decay, the integral J °° xe 2ex \<f>j(x)\ 2 dx con- 
verges, and hence the Hankel operator Tj with symbol e £X (pj(x) is bounded. We decompose 
<pj = + i^scfrj so that we can work with the self-adjoint Hankel operators and Tc Se j )j ; 
so by theorem 2.1 of [44, p. 257], there exist linear systems (-A' j: B^Cj) and (-A",B",C") 
with input and output spaces C, and state space H, and all operators bounded, such that 



! 3ty,-(a;) 



C'e 



-xA' A 



B'j and e ex Q4>j(x) 
(-Aj, Bj,Cj) = 



C'e 



-xA'J 



B'-'; then we introduce the linear system 







" 




B 'i 


(- 





A "\ 


1 


l B "\ 



(5.12) 



which has scattering function e £X (j)j(x) = Cje xA ^Bj. Hence we can introduce the linear 
system 



{-A,B,C) = {- 












el + A 



N 



B 1 


L 







B 



N 



Ci 


L 



c 



N 







(5.13) 



where A : H 2N H 2N ', B : C N -)> H 2N and C : H 2N are bounded linear operators. 

Since ^R(A^,^) h n > e(£,£)h n f° r au £ £ Lemma 2.1 and Proposition 2.2 show that 

(—A, S, C) is a (2, 2) admissible linear system. Evidently (—A, B, C) realises $, and we can 
likewise realise \E» by a (2, 2) admissible linear system. 

□ 

By taking a = 0, 7 to be a negative proper rational function and 1//3 to be a positive 
polynomial on (0, 00), one can produce solutions of (5.2) that satisfy the hypotheses of Theorem 
5.1(ii). 

Remarks 5.2 (i) Whereas Theorem 5.1 does not give an explicit form for the admissible linear 
system (—A, B, C), we can produce one explicitly in several important cases; see section 6 and 
[9, 10, 11]. 
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The next step is to widen the discussion from rational functions on C to meromorphic 
functions on algebraic curves, which involve a spectral parameter. Krichever and Novikov 
introduced the notion of a spectral curve for a family of commuting differential operators [38] . 

Definition (Algebraic family) Let (K, d/dp) be a differential field and let Uj(t,p) be Mn 
matrices with entries in K(ti, . . . , t n ) and let Lj = — Uj(t, p). We write t = (ti, . . . ,t n ) for 
short. 

(i) Say that Lj form a commutative ensemble if [Lj, Lk] = for all j, k. 

(ii) A commutative ensemble is said to be algebraic if there exists a Mjy matrix function 
W(t, p) with entries in K(ti, . . . , t n ) such that [Lj, W] = for all j. In this case the spectral 
curve is 



Lemma 5.3 (i) The curve £ is independent of t and every point [i of £ lies in some finite 
algebraic extension K[/z] of K. For each such fj,, there exists an eigenvector ^(t, p) of W(t, p) 
with entries in K(£, u). 

(ii) Let \I/(i, p) be an eigenvector that corresponds to a simple eigenvalue ofW(t, p). Then 
Lj^(t,p) = A(p,t)*(p,t) for some X(p,t) G K(t,fi). 

(Hi) Suppose further that [A\ is isomorphic to a commutative subalgebra of Mjv(K(t)) 
and that W(t,p) has N distinct eigenvalues as in (ii). Then every [P\ G |_AJ has eigenvalues 
in K(ii, . . .,t n )[p]. 

Proof, (i) We note that d/dtj det(/iijv — W(t, p)) = 0, so det(//7jv — W(t, p) is a polynomial 
in n with coefficients in K, which are integrals of the flows generated by the Lj. Hence each 
(/i,p) on £ is associated with a maximal ideal P of K[/i]; furthermore, K[/i]/P is a finite 
algebraic extension of K by the weak Nullstellensatz theorem [5]. Each eigenvector \& of 
W(t, p) has entries which are rational functions of /i and the entries of W(t, p) hence belong 
to K(tt, . . . ,t n )[fj]. 

(ii) We observe that WLjty = /J,Lj^, hence Lj^ = Xj^ for some Xj G K(t, n). 

(iii) Suppose that #(t,p) satisfies Lj^(t,p) = and *(0,p) = /, hence W(t,p)^(t,p) = 
^(t, p)W(0, p) The columns of the matrix ty(t, p) give eigenvectors of each [P\ G |_AJ, so 
each [fj may be written as a diagonal matrix with respect to the common basis of eigen- 
vectors with entries that are eigenvalues of LPJ- The operators [P\ and W(t,p commute in 
Mn (K(ti , . . . , t n ) and have a common basis of eigenvectors, hence we can write \_P\ \&j (t, fij ) = 
Xjtyj(t, ji) where ^j(t, ji) has entries in K(ti, . . . , t n )[u], hence Xj G K(ti, . . . ,ijv )[/•*]. 



Lemma 5.3 deals with the local behaviour of eigenvectors parametrized by £. Suppose 
that the poles of (5.2) are simple and that the residue matrices are differentiable functions of 
deformation parameters t = (ti, . . . , t n ), so that 




(5.14) 



□ 




(5.15) 
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where trace(t/j) = 0, and consider a family of meromorphic solutions Y = Y(\;t\, . . . ,t n ) of 
the differential equation JdY/dX = Q(\)Y for A complex that also satisfy the conditions of 
Theorem 5.1, and as in (5.2) introduce the kernels 



(1 , ynx+^i^M where Y _ 

( ' x — y 



(5.16) 



Proposition 5.4 Let r(z,t) = det(I + K\J)), suppose that \\K^\\ < 1 for all Siz > xq, and 
suppose that the differential equations 

dY -U 

3 Y (j 1 n) (5.17) 



dtj A — cij 



(i) Then logr(z, t) is given by Proposition 2.4(H) and 



are mutually compatible. 
JL 

dz 



1 ^ £j 

-— logr(z,t) = ^ — logr(z,t) [<Rz>x Q ). (5.18) 
j=i 3 

(ii) Let j be an index such that 5Joj is largest, suppose that 5Joj > 2xo and that 
(JUjY(a,j,t), Y(dj, t)) 7^ 0. Then logr(z, t) has a pole at z = dj/2. 

(Hi) There exists a hyperelliptic curve £ and a complex torus T such that r(A, t) extends 
to £ x T. 

Proof, (i) By a calculation as in [10, Theorem 3.3], we have 

d {JY(x + 2z,t),Y(y + 2z,t)) = / Y(x + 2z,t) Y(y + 2z,t) \ 
dtj x — y \ 3 x + 2z — cij ' y + 2z — cij I 

which decomposes the kernel into a finite sum of integral operators of rank one, and likewise 

1 d (JY(x + 2z,t),Y(y + 2z,t)) = _ A/ Y(x + 2z,t) Y(y + 2z,t) x 

2dz x-y 3 x + 2z-aj ' y + 2z-dj /' l ' ' 

which gives the identity of finite rank operators 



j=i j 

The operator (d/dz)Kr z \ is of finite rank, and hence is trace class if and only if the constituent 
functions belong to {L 2 (0, oo), dx). Now as in Theorem 5.1(h), we choose xo so large that 
/ + is an invertible operator for all z > xq and then compute ^ logr(z, t) = trace((7 + 

/f(zj) _1 Jz-^(z))' so we deduce the stated result. The identity (5.18) asserts that infinitesimally 
translating z is equivalent to the added effect of infinitesimally moving all the tj . 
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In Theorem 5.1, we showed that t(z, t) is given by the Fredholm determinant of a product 
of Hankel operators, and in Proposition 2.4, we expressed ^ log det (7 + 1^ {z) T^ 2 ) in terms 
of the solution of a Gelfand-Levitan equation; thus log t(z, t) is given in terms of the solution 
of a Gelfand-Levitan equation. 

Note that when cij — 2z lies on (0, oo) and Y(a,j) / 0, the function Y(x + 2z)/(x + 2z — a,j) 
does not belong to L 2 ((0, oo); dx), so there is a possible pole for r' (z,t)/r(z,t). 

(ii) We take 2z — dj G C \ (— oo, 0] and compute 

1 d A f°° (JU J Y(x + 2z,t),Y(x + 2z,t)) ^ 
-trace— Ki z \ = — > / dx 

2 dz {z > fW (x + 2z- aj ) 2 



L ^Et^H =0 (7 = 1, ■■-,**) (5.24) 



k=l 

= _ (JU j Y(a j ,t),Y(a j ,t)) +m (j8 _> fl . /2)j (5 . 22) 

so (d/dz)K( z ) has a simple pole at aj/2. By (3.37), (d/dz) logdet(7 + K z ) has a pole at aj/2. 

(iii) Schlesinger [51] observed that the system (5.17) is consistent if and only if the family 
of solutions satisfies an isomonodromy condition with respect to infinitesimal deformation, or 
equivalently that a certain family of differential operators commutes. 

Let V 1 be the space of first order differential operators in time parameters t = (t\, . . . , t n ) 
with coefficients in M 2 (C(A, t)), and let 

L " = w L > = m- + irH- w = i,...»), < 5 - 23 > 

Gamier [25] observed that 

n 

- k=l x - ak - 

hence {Lj;j = 1, . . . , n} gives an algebraic ensemble for the 2x2 matrix 

n 

W(\,t) = m(X,t)'[[(\-a j ) (5.25) 

3=1 

which is a polynomial in A. Since trace(M / ) = 0, we observe that 

det(r ? / 2 - W(X, t)) = rj 2 + det W{\, t) (5.26) 

which is independent of t by (5.24). Hence £ = {p = (A, rj) : rj 2 + det W{\, t) = 0} U {(oo, oo)} 
defines a compact hyperelliptic curve, and the corresponding Jacobi variety T is a complex 
torus of dimension g. Now let K(t) be any differential field that contains the entries of W(t, p). 
Every eigenspace of W(t, p) contains an eigenvector with entries in K(t)[\/det W]; in particular 
on each neighbourhood we can take -0 = column[/, h] and replace it with columnfl, h/f], where 
h/f G K(t) [Vdet W]. Moreover, for each p such that W(t, p) ^ 0, there is a simple eigenvalue 
rj such that the eigenspace {4> : W(t,p)ip = rjip} gives a family of line bundles on £ with 
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additional parameter t, which may be identified with an element of the Picard group of divisors 
on £ modulo the principal divisors by Theorem 3 of [53] . Thus each line bundle has a degree 
n £ Z, and the set of line bundles with n = may be identified with T. 

Hence p = h/f gives an Abelian differential on £, and we can form the integral 

^(p,p ;t) =exp / p(C,t)— (5.27) 

which defines the Baker-Akhiezer function ^(p; po; i); this is multi valued, since the expression 
depends upon a multiplicative factor given by the choice of path from po to p. 

Given a simple eigenvalue r\ of W(t, p), each eigenvector ip is also an eigenvector of Lj for 
all j = 1, . . . , n. Furthermore, by translating A infinitesimally, we make a linear flow in the 
group of divisors, Lq generates rectilinear motion in T. Thus we can extend the tau function 
to 

r(A,i) = det(/ + X ( ( * ) ) ) (t eT,p= (A, rj) G £). (5.28) 

□ 

Remark 5.5 To recover the usual form of Schlesinger's equations [22, 25, 31, 51] one sub- 
stitutes tj = cij after differentiating, and considers the residues at each of the poles. By 
Schlesinger's results, as interpreted in [31], there exists a multi- valued and locally analytic 
complex function rg(ai, . . . , a n ) on 

{(oi, . . . ,a n ) : cij ^ a k ;j,k = 1, . . . ,n} (5.29) 

such that 

d log ts = ^2 t T &ce{UjUk)dlog(aj — cik) (5.30) 

j,k:j<k 

as an identity of differential one forms, so that 

< ** a >- ak 

This contrasts with (5.18), and indicates that translation has a different role for the two 
versions of the tau function. See [14]. 



6. The differential ring associated with the Painleve II equation 

In this section we consider a linear system which is important in random matrix theory. 
Whereas the state ring S is finitely generated, the linear system is not integrable in the sense 
that r does not emerge from C(x) by successive Liouville integrations. Let H(w,v;x) be a 
Hamiltonian which is rational in the canonical variables (w,v) and a meromorphic function 
of time x, and let (w(s), v(s)) be solutions of the canonical equations of motion, and suppose 
momentarily that these are meromorphic functions of s. Then the corresponding tau function 
is 

t(x) = exp / H(w(s),v(s); s) ds, (6-1) 
Jo 
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where the integral is taken along an orbit in phase space; so the value of r is locally independent 
of the path of integration, provided the path avoids poles. 

The Hamiltonians which arise on random matrix theory have additional properties which 
are described in the following result, which is a variant of Theorem 1 in Okamoto's paper [46]. 

Proposition 6.1 Suppose that the Hamiltonian H(w,v;x) is rational in x, a polynomial 
in v, and a quadratic polynomial in w, let u he the potential that corresponds to r. Then 
K = C(x,v,v') is a differential held with respect to d/dx under the canonical equations of 
motion that contains u. 

Proof. We write H = A(v,x)w 2 + B(v,x)w + C(v,x). Then the canonical equations are 
cE = 7Fw anc ^ lie = ~^5v~ Hence C(x)[w, v] is a commutative and Noetherian differential ring 
for the derivative ^ = ^ + 7^7 ~~ l^"§w- Using the special form of the Hamiltonian, we 
have 

~*(dA/v'-B\2 dB(v'-B\ dCi v'-BfdA. dA\ (dB . dB\ 

(6-2) 

so v" = f(x,v,v') where / is rational in x and q and quadratic in v' . So C(x,v)[v'] is a 
differential ring for d/dx. Likewise, the potential that corresponds to r is 

dH dA(v'-B\* dB(v'-B\ dC , . 

hence there exist E,F,G G C(x, v) such that Ev 12 + Fv' + G = u, so K is a differential field 
containing u. 

□ 

Okamoto [46] has shown that each of the Painleve transcendental differential equations 
Pi, ... , Pyi arises from a Hamiltonian as in Proposition 6.1, and r is meromorphic on a suitable 
covering surface. In particular, for x G C and a complex constant a, let 

1/ X \ ^ / X \ / x \ x^ 1 
Hu(w,v;x) = -[w- -J + + -j\w--j -av + —. (6.4) 

Proposition 6.2 Under the canonical equations of motion with Hamiltonian Hu, 

(i) v satisfies Pn : v" = xv + 2v 3 + a and the corresponding r function is 

/ I f°° \ 
r(x)=exp^— - y (5 — x)v(s) 2 dsj ; (6.5) 

(ii) w satishes w'" + 6ww' - (2w + xw') = and U(x, t) = (3t)" 2/3 w(3"H- 1 / 3 x) satishes 

d 3 U 2U dU 1 dU _ 

lW + l^~dx~~9lH~ ^ ' 

Proof, (i) The canonical equations of motion are satisfied in the polynomial ring S a = 
C[x,v,w] with the derivatives 

dv o dw 

— = —w — v and — — = (2w — x)v — a. (6.7) 

dx dx 
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Hence K = C(x, v, w) is a differential field, and by Lemma 7.1 the potential is u = v 2 , which 
belongs to K. We deduce that v satisfies Pjj. 
(ii) Now w satisfies 

„ o a(a + 1) + w' — (w') 2 

K 2 : w" + 2w 2 -xw+ v ; = 0. 6.8 

2w — x 

One can then verify that U satisfies KdV; see [2] for further discussion. 

□ 

Proposition 6.2 indicates that we have potentials linked by the Miura transform, as in 
(2.37), and the next step is to introduce suitable linear systems and operators to solve the 
differential equation Pjj. We introduce Airy's function Ai(x) = e l € x + l £ / 3 d^/(2n), which 
satisfies Ai"(x) = xAi(x). Let (j)(x) = Ai(cc) and let £ = </>'/</>; then S = C[x,4>(x),£(x)} is a 
differential ring with respect to d/dx. In the context of Theorem 6.3(iii) below the integrable 
kernel 

o2/ ^ Ai(x)Ai / (y)-Ai / (x)Ai(y) 

R (x,y) = — — (6.9) 

X y 

is known as Airy's kernel, which is associated with soft edges of eigenvalue distributions [59]. 
The Fredholm determinants of R 2 lead to a solution of the Painleve II nonlinear differential 
equation. Ablowitz and Segur solved Pjj by a slightly different method, Borodin and Deift 
[12] obtained a solution by considering a matrix Riemann-Hilbert problem involving (6.9) and 
we include the proof of Theorem 6.4(iii) to illustrate the general theory of linear systems. 

In previous sections we started from an admissible linear system and produced a Hankel 
integral operator F^. In this section we begin with a technical result which realises a typical 
Hilbert-Schmidt Hankel operator from an explicit linear system (—A, 5, C) chosen for 
(f>. Here A is defined on V(A) = {/ <E L 2 (0,oo);f G L 2 (0,oo)} and C is bounded on V(A). 
Suppose that 4> and ip are continuous functions on R such that / °°(1 +t)(\(j)(t)\ 2 + \i^(t)\ 2 ) dt < 
oo. Then we let H = L 2 (0, oo) and introduce the operators 

A : /Or) -> -f'{x) f E V(A); 
B : ^ (f)(x)p; 
E : (3 i-> 4>(x)f3; 

C:g(x)^g(0) (g G V(A)), (6.10) 

so that 4>{x) = Ce~ xA B and ij){x) = Ce~ xA E. We introduce the operators on H given 
by R x = J™ e- tA BCe~ tA dt and S x = e~ tA ECe~ tA dt. In terms of Lemma 2.1, the 
cogenerator V is unitarily equivalent via the Fourier transform to the coisometry on the Hardy 
space H 2 on the upper half plane 

v . /(,) „ (1 - am - m (/eff2) , (6 . n) 

1 + IZ 

so is the shift. This is consistent with Beurling's canonical model of a linear system in [7]. 
We also introduce the observability Gramian Q x = e~ M+ C^Ce~ tA dt and we observe that 
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Q x is the orthogonal projection Q x : L 2 (0, oo) —> L 2 (x, oo). We consider the Gelfand-Levitan 
integral equation (2.8) where T(x, y) and $(x + y) are 2x2 matrices, and 



$(x) = 



^(x) 





(6.12) 



Lemma 6.3 (i) There exists 5 > such that for all n £ C such that < 8, the operator 
I — jj?R x S x has inverse G x G B(f7) and the matrix function 



f(x,y) = 



^Ce- xA G x S x e-y A B —Ce~ xA G x e~ yA E 
Ce - xA -yA B _ ^Ce- xA R x G x S x e- yA B u.Ce- xA R x G x e- yA E 



(6.13) 



satisfies the Gelfand-Levitan equation (2.8). 
(ii) The determinants satisfy 



det(J - n 2 R x S x ) = det(J - ^Y^Y^] 



and 



1 d 
H dx 



trace T(x,x) = --^ logdet(i" - fJ^Y^T^] 



(6.14) 
(6.15) 



Proof, (i) We introduce 

A = 



'A 


" 


, B = 


' B 


" 


, c = 


" 


c~ 





A 





E 


c 






(6.16) 



and follow the computations of Proposition 2.3(i) to find T. 

(ii) We observe that R x f(z) = <fi(z + u)f(u) du, so R x is a Hilbert-Schmidt operator, 
and Ro is the Hankel operator r^; likewise S x is Hilbert-Schmidt; hence R X S X is trace class. 
The identity (6.14) follows from Proposition 2.3. 

Whereas R X S X is not a differentiable function of x, and we cannot adopt the direct ap- 
proach of Proposition 2.3(iii), we can differentiate the Hankel product 

d r°° 

— (j)(2x + s + u)ip(2x + t + u) du = -2(/)(2x + s)ifj(2x + t), (6.17) 
dx Jq 



so 



dx^ 4 '^^^ ~ 26 



-2xA 



BCe 



-2xA 



So, 



(6.18) 



where the right-hand side is a rank one and bounded linear operator. We recall from [48] the 
following identities regarding the shift and Hankel operators 



6 x ^ g 



= Qx, e- xA e- xAf = I, R e- xAi = e~ xA R Q (6.19) 



and the following special identities which may be checked by looking at the kernels 



Rx — RoQx, e xA Rq — R x e 



-xA j 



_ — xA p _ — X 

1 4>(x) - e R e 



-xA^ 



(6.20) 
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Hence we can differentiate using (6.15), obtaining 

— logdet^-^r^r^) 



= 2 M 2 trace((/ - ^T^T^) 1 e - 2xA BCe- 2xA S ) 



= 2/j 2 Ce- 2xA S {l - v 2 e- xA R e- xA e~ xA S { 
= 2^ 2 Ce- xA S x e- xAi (I - ^ 



oe 



?A^\- 1 -2xA 



D 



2-xA 



Roe 



-iA T — xA 



S e- xAi ) 



A i \- 1 ~2xA 



B. 



(6.21) 



We now use the identity K(I + LK) 1 = (/ + KL) 1 K to shuffle terms around, and obtain 



= 2^ 2 Ce- xA S x (I - v 2 e- xAi e- xA R e- xAi e~ xA 



Soe-^r'e-^e-^t 



V B 



= 2^ 2 Ce- xA S x (I - ^QxRoQxSoe-^^Qxe-^B 
= 2/j 2 Ce- xA S x (I - yPRxSxY^e-^B 
= 2^ 2 Ce- xA (I - ^SM^S^-^B; 

which is a multiple of the top left entry of T(x, x), and likewise 

2ii 2 Ce- xA R x {l - fi 2 S x R x )- 1 e- xA E = 2 l x 2 Ce- xA R x G x e- yA E; 

as in the bottom left entry of T(x, x) so we obtain the expected result 

d 



dx 



logdet(7 — /j, R x ) = / utraceT(x, x). 



(6.22) 



(6.23) 



(6.24) 



□ 



Let {—A, B, C) be as in Lemma 6.3, where T(x, y) and $(x + y) are 2x2 matrices, and 



W(x,y) V(x,y) 
V(x,y) W(x,y) 



$(x) = 



T(x,y) = 

where the entries of T are given by 

W(x, y) = iiiCe- xA {I + n 2 R 2 x y l R x e~ vA B 

and 

V(x, y) = -Ce- xA (I + v?RlY x e- yA B. 



o 4>{x) 

<f>(x) 



(6.25) 



(6.26) 



(6.27) 



Theorem 6.4 (i) There exists S > such that for all \i G C such that < 8, the operator 
I + /J. 2 R 2 is invertible, and T satisfies 

/>oo 

T(x,y) + $(x + y) + / T(x, z)$(z + y)dz = 0. (6.28) 

J X 

(ii) Suppose further that (—A, B, C) is (2, 2) admissible. Then the determinant satisfies 

(6.29) 



V(,,x) = J-iLlog drt(/ + ^» ) 



2%\xdx det / — i/jT^^A ' 
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W(x,x) = ^^et(I + u 2 ^), (6.30) 
1 d 



2ifi dx 

hence 

d -W(x,x) = -V{x,x) 2 . (6.31) 



2i/j, dx 

(in) In particular, let <p(x) = Ai(x/2). Then v(x) = V(x, x) satisfies Painleve's equation 

d 2 

P n -j-ju(i) = xv(x) - 8u. 2 v(xf (6.32) 

and v(x) x — Ai(cc) as x — >• oo. 

(iv) The entries of T(x, x) all lie in the differential ring C[x, v, v'], and the potential is 

u(x) = -8n 2 v{x) 2 . (6.33) 

(v) The cumulative distribution function of the Tracy-Widom distribution [56] satisfies 

F 2 (x) = det((/-r 2 ^ ) /4). (6.34) 

Proof.(i) This follows from Proposition 2.3. 

(ii) The first part follows from Proposition 6.2(h), while the second part is as in Theorem 

2.6. 

(hi) First, note that V(x,x) = -Ce~ xA (I + fi 2 Rl)- 1 e- xA B where Ai(x/2) = Ce~ xA B, 
so V(x, x) is asymptotic to — Ai(cc) as x — >■ oo. 

It follows from the Gelfand-Levitan equation that 

/»oo /»oo 

V(x,y) + (f>{x + y)+n 2 / V(x, z)<f>(z + a)<f>{a + y) dzds = 0. (6.35) 

J X J X 

Let L = (-^ + j^) 2 - 2±a an d 0( x ) = Ai(x/2), so that L(f)(x + y) = 0. Also from Airy's 
equation, we obtain 

/>oo 

y - z 



(j)(z + t)(f)(t + y)dt = 4((f) / (z + x)(f)(x + y)-(j)(z + x)(t)'(x + y)), (6.36) 
and by repeatedly integrating by parts, we can reduce (6.27) to the expression 

^ />oo 

LV(x,y)-4/i 2 ( — J V(x,z)(j)(z,x)dz)(j)(x + y) 

+ n 2 I I LV(x,z)<f>(z + s)<f)(s + y) dzds = 0. (6.37) 



Now we have 



/>oo 

W(x,y) + i(i V(x,z)(f>{z + y)dz = 0, (6.38) 
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so 



— / V(x, z)6(z + x)dz = j-W(x,x) 

ax J x ifx ax 



dx . 

= 2V(x,x) 2 . (6.39) 
On multiplying (6.38) by — 8fj?V(x, x) 2 and using uniqueness, we deduce that 

LV(x, y) = -8/j, 2 V(x, x) 2 V{x, y), (6.40) 

and on the diagonal we have 

d 2 

P n -^V(x,x)-xV(x,x) = -8/j 2 V(x,x) 3 . (6.41) 

(iv) From (6.32), we see that C[x, v, v'\ is a differential ring. Note that w'(x) = —2i^v(x) 2 
so the differential equation gives 

w(x) = 2ifi(-xv(x) 2 + v'{x) 2 + 4^ 2 v(x) 4 ), (6.42) 

which are all elements of C[x, v, v'], hence the entries of T(x, x) are all in C[x, v, v']. 

(v) With fj, = i/2, the potential gives rise to the Tracy Widom distribution function 

F 2 (x) = exp(-2- 1 f (s-x)u(s)ds) (6.43) 

J X 

that is associated with the soft spectral edge of the Gaussian unitary ensemble; see [55, 56 
(1.17)]. 

□ 

The results of this section extend to differential equations which arise in unitary matrix 
models. For each £ = 1,2,..., there exists a real solution Ag(x) to the generalized Airy 
equation 

Af £ \x) + (-lYxA e (x) = (6.44) 

such that Ai(x) — > exponentially as x — >■ and Ag(x) is algebraically damped and oscillating 
as x — >■ — oo. In [16], Crukovic, Douglas and Moore select a suitable curve C in the integral 

1 f ( w 2e+1 \ 
A e (x) = - / exp i [wx + — — — )dx (6.45) 

7T Jq V At + 1/ 

and use the stationary phase method to prove the existence and uniqueness of (f>. 
Definition (i) Let 4>{x) = Ag(x/2), and let the generalized Airy kernel be 

t - Z 



1 <b (2e - k) (tU {k - l) (z) - 6 {k - l) (tU {2e - k) (z) 

K (x) (t, z) = 2 2 ^ ^(-i)fc-i-^W {Z) J^ {Z) 



(6.46) 



k=i 
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which reduces in the case £ = 1 and x = to the standard Airy kernel as in [59]. 

(ii) Let J be a finite union of open sub intervals of (0, oo) and let Pj be the orthogonal 
projection on L 2 (0, oo) given by multiplying functions by the indicator function of J. A point 
configuration on J is a function v : J — > {0, 1, 2, . . .} such that {x G S : u(x) > 0} is finite 
for all compact subsets S of J; write u(S) = Ylxes v(x). Let Conf(J) be the set of all point 
configurations on J with the cr-algebra generated by the functions u, so v \-t u(S) is a random 
variable on Conf(J). We regard u(S) as the number of random points in S. 

(iii) For the remainder of this section, we temporarily change notation, so that we conform 
with the standard references. Gelfand and Dikii [26] introduced a sequence of polynomials 
Rj [u] G R[ix, u', u", . . .] by the recurrence relation 

R' j+1 = 4 _1 i2j" - uRj - 2~ 1 u'R j (6.47) 

and the sequence begins with 

Rq = 2 -1 , R 2 = -4 _1 ix, R 2 = W^^u 3 - u"). (6.48) 



Proposition 6.5 (i) Then = , so defines a trace class operator on L 2 (0,oo) 
such that < Kt~-\ and K<~\ — >■ as x — >• oo. 



L (a .) ana /v (a .) 

(iij Let w be the differentiated quotient of determinants 



then Jet u(x) = i/(x) + i>(x) 2 , so that it is the Miura transform of v. Then v satisfies the string 
equation 

p n ■ " vR t [u] + (-l) e+1 2- 2e xv(x) = 0, (6.50) 

which reduces in the case I = 1 to Pn . 

(iii) Then there exists a probability measure on Conf((0, oo)) such that 

(—l) n d n 

E(n;J) = { -^~^ x= det(I - XK {xo) Pj) (6.51) 
equals the probability of {u : v(J) = n} for all J and all n = 0,1, ... . 

Proof, (i) Since is real and rapidly decreasing, the Hankel operator r^ (a:) is self-adjoint and 
Hilbert-Schmidt, hence L 2 ( } > is trace class. To prove that = F 2 ( , we multiply the 
right-hand side of the claimed identity by (t — z) and use the differential equation to obtain 

{-iy + \t- Z ) 



<j>(t+s)<j)(s+z)ds= / (<j)^ ) {t+s)^{s+z)-^{t+s)^ 2i \s+z))ds. (6.52) 



>2x J2x 

After integrating by parts I times, we obtain the formula 



poo 

K {x) (t, z)= 4>{t + s)(f)(s + z)ds. (6.53) 

J2x 
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(ii) This follows from Theorem 2.6. Zakharov and Shabat [62] show that the Gelfand- 
Levitan equation leads directly to a system of differential equations, which in this case gives 
the generalized Painleve equations From the recurrence relation, we obtain the sequence of 
differential equations that starts with 

Pff : v" + 2v 3 -xv = 0, 

P\f : v [A) - 10v"v 2 - \Q{v') 2 v + 6v 5 + xv = 0. (6.54) 

(iii) For all xq > sufficiently large, the kernel satisfies < Kr x \ < I for all x > xq 
and by Mercer's theorem PjKPj is a trace class kernel for any compact subset J of (0, oo) 
as in [55]. Hence Ki x \ defines a determinantal random point field on (xo,oo) by Soshnikov's 
theorem [Theorem 3, 53]. We refer the reader to this paper for more precise definitions of 
the construction of the measure such that E(n, J) is the probability of {v G Conf((x , oo)) : 
v{J) = n}. 

□ 

The string equation (6.50) arises from a double scaling of a probabilistic model of unitary 
random matrices, so one would expect Ki x \ to have a probabilistic interpretation in the style 
of Tracy and Widom. The system of differential equations (6.54) is associated with scaling 
solutions of the Zakharov-Shabat hierarchy, and is sometimes known as the Painleve II hierar- 
chy; although there are competing notions of this concept. One can also introduce a hierarchy 
of differential equations from the system from the dispersive water wave equations, the first 
terms of which reduces to the system of linear equations which that Jimbo and Miwa [30] use 
to linearize P//; but at the time of writing it is not clear as to whether the various approaches 
lead to equivalent hierarchies. 

7. The differential ring of a periodic linear system 

In this section we obtain analogues of Theorem 4.2 for periodic groups. For periodic 
and meromorphic u, the differential equation — i/j" + uip = is known as the complex Hill's 
equation. In subsequent sections, we consider special periodic linear systems such that u is a 
function of rational character on the cylinder or u is doubly periodic and of rational character 
on some elliptic curve T . 

For periodic linear systems, the defining integral for R x in Lemma 2.2 does not converge, 
and the contour integral for Rq in Lemma 4.1 is inapplicable; nevertheless, we can adapt a 
result of Bhatia, Dacis and Mcintosh discussed in [8] and otherwise construct R x satisfying 
(1.10). 

Lemma 7.1 Let B be a trace class operator and C he a hounded operator on H, and let 
(e~ M )teR he a bounded Cq group of operators on H such that the spectrum of A does 
not intersect the spectrum of —A. Then there exists a solution to the Lyapunov equation 
— -^R x = AR X + R X A such that AR + R A = BC and R x is trace class for all x G R. 

Proof. The main problem is to find E such that EA + AE = BC. By a theorem of Sz-Nagy, 
the group (e~ tA ) is similar to a group of unitaries, so there exists an invertible operator S and 
a unitary group (U t )ten such that e~ tA = SUtS^ 1 . Hence the spectrum of the skew adjoint 
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operator A lies on iH and is a closed subset. By hypothesis, there exists 8 > such that the 
spectra of A and —A are separated by 8 and <j(A) U a (—A) does not intersect (—8,8). By 
Plancherel's theorem, we can construct an integrable function / such that /(£) = l/£ for all 
£ G R such that |£| > <5. Then the integral 

/oo 
e- xA BCe~ xA f(x)dx (7.1) 
-oo 

has a weakly continuous integrand in the trace class operators, and is absolutely convergent 
with 

/oo 
||S|| c i||C|| B(H )M 2 |/(x)|da; (7.2) 
-oo 

hence E is trace class. Using the spectral representation of U t , one can show that AE + EA = 
BC . Next we introduce R x = e~ xA Ee~ xA which gives a one parameter family of trace class 
operators such that —^j^ = AR X + R X A. On verifies this identity on T>(A) and then observes 
that both sides are trace class. 

□ 

Definition (Periodic linear system) Let (e~ xA ) x<E n be a uniformly continuous group of opera- 
tors on H such that e 27TA = I and A is invertible. Suppose further that B and E are trace class 
operators on H, and that C is a bounded linear operator on H, such that AE + EA = BC. 
Then Sqo = (—A, B,C; E) is a periodic linear system with input, output and state spaces all 
equal to H. Whenever we define a parametrized family Et of periodic linear systems, the input, 
output and state spaces are taken to be fixed; furthermore, A is taken fixed in the family. 

We let C = C/ttZ be the complex cylinder formed by identifying w~zifz-wG 7rZ; 
we can choose equivalence class representatives in the strip {z : — 7r/2 < !Kz < 7r/2}; then we 
identify each 7r-periodic / : C — >■ X with a function / : C — > X. Let Cc = C[sin2z, cos2z] and 
let Kg = C(sin2z, cos2z) be the field of trigonometric functions, which consists of functions 
of rational character on C in the sense that the elements are rational functions of t = tan z. 
The space of entire tt periodic functions on C may be identified with the space of holomorphic 
functions He on C, which is differential subring of the meromorphic functions Mc on C. 

Definition (Tau functions and operators) Adjusting the definitions of section 3 in a natural 
way, we let <&(x) = Ce~ xA B be the operator scattering function so that <f>{x) = trace $(x) 
is the (scalar) scattering function and let R x = e~ xA Ee~ xA , then we introduce F x = (I + 
e~ xA Ee~ xA )~ 1 , and t^x) = — detF^, then let u(x) = —2-^logT 00 (x) be the potential. Let 
Spec(A) be the spectrum of A as an operator, and introduce the periodic linear system 

S A = (-A, (XI+A)(XI-A)- 1 B, C; (\I+A)(\I—A)~ 1 E) (A e (Cu{oo})\Spec (A)) (7.3) 

and its accompanying tau function t\. Whereas the linear system has infinitely many inputs 
and outputs, we impose severe spectral conditions on A in lieu. We also introduce the (non 
commutative) algebra S = Kc{I, A, BC, F x }, and then let A be the subring of S spanned by 
A ni and by the ordered products A ni FA n2 . . . FA Ur for rij G N. We also introduce [_ • J : S — >■ 
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Mefc 1 ) : [P\ = Ce- xA FPFe- xA B. Let [AJ = {[PJ : P G A} and A = {traceLPj : P G A}. 
Note that A is analogous to the differential ring generated by the potential u. 

Theorem 7.2 Let (— A B, C; E) be a periodic linear system. 

(i) Then t\(x) is holomorphic except at fixed singularities; so x ^ t\(x) is entire, while 
A ' V T\(x) is holomorphic on C \ Spec(A); 

(ii) 4>{2x) G Cc, and S is a complex differential ring for (—A, £?, C; E) and for E^; 
(Hi) |_AJ is a complex differential ring on C; 

(iv) the derivatives of the potential belong to Mc and to Aq. 

(v) Ife- A7T / 2 Ee- AlT / 2 = -E then T(x,y) = -Ce~ xA F x e- yA B satisfies 

^ fX + TT/2 

$(x + y) + T(x, y) + - J T(x, z)<5>{z + y) dz = 0. (7.4) 

Proof, (i) First we show that A is an algebraic operator. By periodicity, the group (e~ xA ) xe n 
is bounded and hence by Sz-Nagy's theorem, e xA is similar to a unitary group on H, so A is 
similar to a skew symmetric operator. By uniform continuity, A is bounded, and hence has 
spectrum contained in {— iN, . . . , iN} for some integer iV; see [20]. Consequently, there exists 
a monic polynomial p such that p(A) = 0. 

Hence A is an invertible algebraic operator, so as in (4.5), A^ 1 is a polynomial in A and 
(AJ + ^XA/-^)- 1 G S for all A in the resolvent set of A Observe that (XI + A)(\I - A)~ l is a 
polynomial in A with coefficients that are rational functions of A, and holomorphic except when 
A is in the spectrum of A; in particular it is holomorphic on {A : |A| < 1} U {A : |A| > ||A||}. 
Hence t\ is a holomorphic function of A, except at A G Spec(y4), which is a finite set. 

(ii) We also introduce polynomials pj for each point in the spectrum of A such that 
Pj(ik) = 8jk for k = —N, . . . , N, and since A is similar to a skew operator, we deduce that 

N 

e ' xA = E PM)z- ijx , (7-5) 

so <&(x) = Ce~ xA B is a trigonometric polynomial with coefficients in c 1 and of degree less 
than or equal to N. Hence 4>(2x) is 7r-periodic. 

By (7.5) and (7.1), the operator E belongs to S and hence R x = e~ xA Ee~ xA also belongs 
to S. Hence we have 

—R x = - e -**AEe- xA - e~ xA EAe~ xA = - e - xA BCe' xA 
dx 

and so AF + FA- 2FAF = F e - xA BCe- xA F, hence 

dF 

— = AF + FA- 2FAF- (7.6) 

dx 

so S is a differential ring for (—A, B, C). 
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(ii) From (7.6), we have the product rule 

[P\ [Q\ = [P{AF + FA- 2FAF)Q\ , (7.7) 

and just as in Theorem 3.3 

A [ P \ = [a(J - 2F)P + ^ + P(/ _ 2F)A\ . (7.8) 

As in Lemma 3.2, 

LAJ = span c {Ce- xA FA ni Fe- xA B,Ce- xA FA ni FA n2 . . . FA n 'Fe~ xA B; rij G n} (7.9) 
is a differential ring. 

(iii) Since e~ xA is an entire operator function, we deduce that is entire, and tt periodic 
since Too(x) = det(I + e 2xA E) and e 2nA = I. When t^x) ^ 0, we have 

4- logdet(/ + e- xA Ee~ xA ) = -trace ((/ + e-^Ee'^^e-^iAE + ^AJe"^) 
ax 

= -trace((7 + e-^Se-^J-^-^SCe - ^) 
= -trace (Ce-^^ (7 + e- xA Ee- xA Y^e- xA B) 
= -trace(Ce- xA Fe- xA B), (7.10) 

and hence 

r/ 2 

u = -2— log det(/ + e-^Se"^) 
dx l 

= -4tvace Ce- xA FAFe~ xA B 

= -4traceLAJ; (7.11) 

so u belongs to A = {trace \P\ : P G A}. Likewise, the derivatives belong to A since 
|_AJ is a differential ring. 

(iv) One can verify this by direct computation, and the crucial identity is 



e~ zA BCe~ zA dz = [-e- zA Ee~ zA ] * +7r/ = 2e~ xA Ee~ xA . (7.12) 



□ 



Remark If (7r/4)||$|| 0O < 1 in Theorem 7.2(v), then 



' )2 ;T(x,y) - -f^T(x,y) = -2(^-T(x,x))T(x,y), (7.13) 



dx 2-y^y> op y > y '~ ^dx 

as one can prove by substituting in the integral equation. This motivates the definition of u 
as the scalar potential, since u(x) = — 2^traceT(x, x) by (7.10). 

If we assume more commutativity, the proofs simplify and the results become stronger. 
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Corollary 7.3 Suppose further that ABC = BCA, and let E = 2- 1 A- l BC. 

(i) Then R x satisfies (1.10) and (1.11); 

(ii) (—A, B, C) is finitely generated, since the algebra S is commutative and Noetherian, 
and is complex state ring for (—A, B, C) on C = C/27rZ. 

(in) Let p he a maximal ideal in S . Then each y G S/p may be identified with a rational 
function on a compact Riemann surface y. 

Proof, (i) Since A^ 1 and C are bounded and B is trace class, E is also trace class. Now 
R x = e~ xA Ee~ xA is an entire and trace class valued function, and using commutativity, one 
checks that Lyapunov's equation (1.8) holds. Unlike in Lemmas 2.1 and 6.1, we do not assert 
that the solution is unique. 

(ii) Here e~ xA is a polynomial in A, e tx and e~ lx , hence e~ xA and likewise R x belong to 
K C ' [/, E, A]. Observe that the set S = {(I + e - xA Ee~ xA ) n : n = 0, 1, . . .} is multiplicatively 
closed and does not contain since I -\- e~ xA Ee~ xA is invertible in the Calkin algebra ofB(H) 
modulo the compact operators on H. Hence we can identify S with the ring of fractions of 
Kc [A, BC] modulo S. There is a natural surjective ring homomorphism Kc [Xi, X 2 , X 3 ] — > S 
given by X\ \-t A, X 2 i-> BC, X 3 — >■ F x , so by Hilbert's basis theorem, S is Noetherian as a 
commutative ring. 

(iii) An ideal p of S is maximal, if and only if {p} is closed in the prime spectrum Spec(S), 
with the Zariski topology. In this case the field S/p is isomorphic to a finite algebraic extension 
of Kc , by the weak form of Nullstellensatz as in [5] . 

Now we produce a finite algebraic cover of P 1 . For each a G S/p there exist Oj G Kc, 
with a n 7^ 0, such that X^=o a ? Qj = 0- By changing variables to t = tanx/2 we replace 
trigonometric functions by rational functions, and multiplying by a suitable polynomial in t, 
we can introduce qj(z) G C[z] such that X)j=o 1j(^) aJ = 0> thus (a,t) is associated with the 
curve {(w,t) : Yl^olji^)^ = 0}j which determines a Riemann surface y which covers P 1 
finitely. 

Each non-zero /i in the spectrum of the compact operator BC is associated with a finite- 
dimensional eigenspace, which is invariant under the operation of A, so we can introduce a 
common eigenvector v of A and BC. Then for each e tx G C there is a maximal ideal P and 
an isomorphism S/p = C given by A \-> ik, BC fj, and R t-> e ~ 2lkx \x /2ik so that p gives a 
point on 

y = {(£, A) : det(/ + \R X ) = 0; e ix = (1 - £ 2 + 2i£)/(l + ^ 2 )}. (7.14) 

□ 

We now consider the tau functions of periodic linear systems (— A, B, C;D). By taking 
traces or forming determinants, we carry out limiting processes which generally take us from 
Kc to Me; the determinant of an infinite matrix of trigonometric polynomials can produce 
rather complicated functions. The scattering function conveys information about the spectrum 
of A, while the zeros of determine the poles of u. This is made precise in the following 
result. 

Proposition 7.4 Let {—A, B, C; E) be a periodic linear system as in Theorem 7.2, and let t\ 
be the tau function ofY,\. 
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(i) G He satisfies log + log + |too(z)| < 2N\z\ + c\ for some c\ > and all z, where N 
is the spectral radius of A. 

(ii) Let (t\) = {z G C : t\(z) = 0} for all A G (—00, 00) U {±00}, which is either empty or 
countably infinite. Every zero of t\ gives rise to a double pole of u\ = —2(logT\)" . 

(Hi) If E has finite rank, then Too is of exponential type and in Cc ■ Conversely, if Too is 
of exponential type, then there exist ctj G C, a G Z and (3 G C such that 

m 

Too(z) = e 2taz+ ? Yl sm2(z - a,) (7.15) 
i=i 

and 

m _ 
^sin 2(z-ctj) 

Proof, (i) Let aj(V) be the approximation numbers of a compact operator V. Then we have 
a n (e~ zA Ee~ zA ) < ||e _i: ' 4 || 2 a n (£^) and hence by a standard bound on the determinant 

00 

\og\det(I + e- gA Ee- zA )\ < c e 2 ^l J>,(£). (7.17) 

i=i 

(ii) If t\(z) = 0, then t x (z + fcvr) = for all k G Z. 

(iii) There exists a projection P of finite rank p such that PEP = E and hence r ao (z) = 
det(7 + PEPe- 2zA P), where Pe~ 22A P is a finite matrix with entries that are in Cc; in 
particular, the entries are functions of exponential growth. Hence from the expansion of this 
determinant, we deduce that there exist c\,ci > such that |too(z)| < cie 2pAr ' 2 ' +C2 for all z. 

Suppose conversely that r is of exponential type, and recall a standard argument of 
function theory. By Jensen's formula, the number of zeros of Too inside a circle of radius r 
grows like c?^ + C4 for some 03,04 > 0, and since Too is also 7r-periodic, we deduce that there 
exists m < 00 such that the only zeros of Too in {z : — vr/2 < ?ftz < n/2} are a\, . . . , a m ; there 
there exists an entire function g such that 

m 

T 00 {z) = e 9 ^Y[sm2(z-a j ), (7.18) 

j=i 

where g is an entire function such that g(z + n) — g(z) = 2-Kit for some t G Z. Since | sin(x + 
iy)\ — >■ 00 as y —> 00, we deduce that \g(z)\ < c^\z\ + ce for some c^,ce > 0, and we finally 
obtain g(z) = 2%az + j3 where a G Z. 

By computing « = — 2(logT 00 ) // , we obtain a potential as in (9.10), which is a rational 
function of e lx and e~ tx . In particular, when m = 1 we have u(z) = 8/ sin 2 2(z — a±), so we 
can rescale this to the familiar case of Csech 2 z for some C. 

□ 

Remark 7.5 The potential (7.16) can be interpreted in terms of a simple model in electrody- 
namics, considered by Sutherland [57]. Consider m fixed unit charges placed at points e lolj on 
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a circular ring, and a further unit charge which has variable position e lx on the ring. Then the 
electrostatic energy of the moving charge is u. In section 10, we show how this can otherwise 
be realised as a limiting case of periodic linear systems with elliptic potentials. 

We now introduce the analogue for periodic linear systems of the Miura transform. Let 
(—A, B, C; E) be a periodic linear system, and let 



T(x,y) = 



W(x,y) V(x,y) 
V(x,y) W(x,y) 



$(x) = 



<j>{x) 
<f>(x) 



(7.19) 



where 



V(x, y) = -Ce- xA (I - {e- xA Ee- xA f)- l e- yA B, (7.20) 

W(x, y) = Ce- xA (I - {e- xA Ee- xA fY 1 e- xA Ee- xA e- yA B (7.21) 

and <f>{x) = Ce- xA B. Then we let r (x) = det(J + e~ xA Ee~ xA ) and n(x) = t (x + vr/2). 
Corollary 7.6 Suppose that e -*A/2 Ee -nA/2 = _ E Then 

(i) T satisfies the Gelfand-Levitan equation 

T(x,y) + $(x + y) + - / T(x, z)$(z + y) dz = 0. (7.22) 

(ii) The entries ofT(x,x) satisfy v{x) = tracey(x,v) and w(x) = trace W(x,x) where 

(7.23) 



id n(x) 

2 dx To \x) 



w{x) = ^ — logr (x)ri(x). 



(iii)Ifro has all its zeros simple and To(zo) = implies to(zo+tt/2) / 0, then v(x) 2 +w'(x) 
has poles of order less than or equal to one. 

Proof, (i) This follows from Theorem 7.2(iv). We introduce the periodic linear system 



-A 
—A 





B 










C 




1 





B 


5 


C 








E 
E -0 



(7.24) 



and follow the proof of Proposition 2.5. 

(ii) This follows from a similar calculation to Theorem 2.6. 

(iii) This is a purely local calculation. The poles of v and w are all simple and occur at 
the zeros of tq(z) and the zeros of tq(z + 7r/2), with no cancellation between these; however, 
the second order poles of v 2 cancel all the second order poles of w'. 

□ 

8. Differential rings related to the KdV hierarchy 

In this section we consider a differential ring S associated with a linear system (— A,B,C), 
and obtain solutions of the Korteweg- de Vries equation and its higher order counterparts. 
First we introduce an infinite collection of time variables tj and consider solutions which are 
parametrized by the infinite torus. 
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We introduce the infinite real torus 

T°° = G R°° : limsup \t 23 -i\ 1/j 0, j -> 00} (8.1) 

which forms a vector space under coordinate wise addition and multiplication by scalars. We 
regard the space variable x as a time parameter t-\. By means of the map (£2,7-1)^1 
YlJLi ^2j—\z^~ , we embed T°° as a linear subspace of the space He of odd entire functions 
that are real on the real axis, where He has the topology of uniform convergence on compact 
subsets of C. The linear functionals £ 2 j-i(f) = /^ 2j ' -1 H0)/(2j — 1)! are continuous, and any 
linear map a : C 9 — >■ span{£j : j = 1,3,...} is continuous and has a continuous transpose 
a* : He — >■ C 9 . The matrix that represents a*, with respect to (z 2 ^~ 1 )'jL 1 and the standard 
basis of C 9 , has finitely many non-zero entries. 

As in [31], for z G C with \z\ > 1 we introduce [^J = {2z~ 2j+1 /{2j-l))f j _ 1 , sot ^ 
is the Sato's shift on odd coordinates. 

For v,to G T°°, the equation t = vs + to gives rectilinear motion. We regard £2.7-1 as 
deformation parameters, and T°° as a subgroup of the formal Lie group {exp Y^JLi hj-id 2 j-i '■ 
hj-i G R} with d 2 j-i = d/dt 2 j-i. 

Now we construct representations of T°°. For a linear system S = (—A, B, C), with A as 
in the following result, we write (i • A) = Y%Li hj-tA 2 ^ 1 for t = [t 2j ^)f =1 G T°°. 

Lemma 8.1 Suppose that either 

(i) A G B(H), or 

(ii) (e~ xA ) xe -R is a bounded Cq group of operators on H. 

Then U(t) = exp(t • A) for t G T°° defines a family of bounded linear operators on H such 
that U(t + s) = U(t)U{s) and 

U(lt\) = (CI + A)((I- A)- 1 (3*C > x ). (8.2) 

As t undergoes a rectilinear motion with velocity v, the resolvent operator satishes 

= (v ■ A)R + R(v ■ A). (8.3) 

ov 

Proof, (i) When A G B(H), the series Y^JLi hj-iA 2 ^ 1 converges, and the addition rule is 
clear. The elementary identity, 

(CI + A)( C I - A)-- = ex P (g (2j + 1)c2j+1 ) , (8.4) 

where the series Y^jLi A 2 i +1 / (2j + l)C 2 - ?+1 converges for \(\ > ||A||, so we can use the left-hand 
side as a definition of the right-hand side for all ( outside the spectrum of A. 

(ii) By a theorem of Sz Nagy, (e~ xA ) xe -R is similar to a unitary group on H, so A is 
similar to a skew symmetric operator A. Let a be the spectrum of %A and p the scalar spectral 
measure of A, and a measurable family of Hilbert spaces so that H = j a H^p(d^) is the 
direct integral decomposition given by the spectral theorem; and for t = (t2j-i)^Li e T, 
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let U(t) = exp(^°^ =1 t2j-iA 2 ^~ 1 ) be the bounded linear operator that operates on as 
multiplication by exp(^°^ =1 t2j-\{— i£) ); then U(t) is similar to a unitary and U(t)U(s) = 
U(t + s) for all s,t£T. One can obtain the identity (7.1) from the spectral theorem. 

In either case, we note that U(t)R x U(t) gives the resolvent operator for E(t), where R x 
is the resolvent operator for E. The differential equation (7.2) follows directly. 

□ 

Definition Here U(t) is the odd deformation group of the linear system S(0) = (—A, B , Cq) 
which is associated with the odd deformation parameters t^j-x- We introduce the family of 
linear systems 

E c (i) = (-A, U([(\)U(t)B, CUit)) (t G T°°). (8.5) 

We allow C : H ->• C and B : C ->• # to evolve through time so that C = C U(t) 
and B = U (t)Bo for some initial Co : — >■ C and I?o : C — >■ H, and correspondingly 
R(x,t) = U(t)R x lI(t). The formulas involving C, B and i? are symmetrical with respect to 
time evolution, since B and C both evolve under the same group. In contrast to Corollary 
7.3, we do not assume that A commutes with BC; that BC here will typically have rank one, 
whereas A will have infinite rank. The operation of -^j- on det(7 — R 2 .) is described by the 
Lyapunov equation (1.10) in the form of the following commutator identity 

w,- < 8 - 6 > 

which is analogous to (19) in [38]. 

Proposition 8.2 Suppose that A is bounded and let F x = (I + i?) -1 . Then 
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—R 




= -2 
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A = span c | A ni , A ni F x A n2 . . . F x A n ' r : m , n 2 , . 

is a differential subring of C°°((0, oo) 2 ; B(H)), and the map |_ . J 

[P\ = Ce- xA F x PF x e- xA B 



. ,n r e n| 

A -)> C°°((0,oo) 2 ;C) 



(8.7) 



has range |_AJ , where |_AJ is a differential ring with pointwise multiplication and derivatives 
d/dx and d/dtx- 

Proof. As in Lemma 3.2, the basic relations are 



A(I - 2F X )P +^P + P(I ~ 2F X )A 



^-[P\= A 3 (I-2F x )P + ^-P + P(I-2F x )A 3 

OT3 L C/C3 

[P\ [Q\ = [P(AF X + F X A - 2F X AF X )Q\ . 



(8.9) 



3.10) 



Indeed it follows from the Lyapunov equation (1.10) that F x satisfies the differential 
equations 



dx 



= AF X + F X A 



2F X AF X , 



(8.11) 
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BF 

= A 3 F X + F X A 3 - 2F X A 3 F X (8.12) 

013 

and hence the derivatives from the first and last factors in (8.8) satisfy 

^-Ce- xA F x = Ce- xA F x A(I - 2F X ), ^-F x e~ xA B = (I - 2F x )AF x e~ xA B; (8.13) 

—Ce- xA F x = Ce' xA F x A 3 (I - 2F X ); ——F x e~ xA B = (I - 2F x )A 3 F x e- xA B. (8.14) 
By applying Leibniz's rule, we deduce that |_AJ is closed under d/dx and d/dt^. Furthermore 

F x e- xA BCe~ xA F x = AF X + F X A - 2F X AF X , (8.15) 

so |_AJ is closed under multiplication, and the product rule (8.10) holds. 

Lemma 8.3 Suppose that CqA 4 : H — >■ C and A a Bq : C — >■ H are bounded. 

(i) Then the scattering function 4>{x;tz) = Coe~ 2tsA ~ xA B satisfies the linearized Ko- 
rteweg -de Vries equation 

(ii) Let v(x, £3) he as in (2.?), so that 

v(x,t) = -C e- xA - t3A3 (I + R)- l e- xA - t3A3 B Q ; (8.17) 
and let u(x,t 3 ) = -2§f . Then 

d 2 

M ( X) t) = -2— logdet(/ + i?) (8.18) 



belongs to |_AJ and satisfies the KdV equation 



^ du d 3 u ^ du 
dts dx 3 dx 



Proof, (i) This follows from a simple computation, 
(ii) We have the following table of derivatives 



dv 
dx 
d 2 v 
dx 2 
d 3 v 
dx 3 

We shall prove that 



2|AI; 

4LA 2 J -8LAi^4J; 

8[A 3 \ - 24[A 2 F X A + AF X A 2 \ + 48[AF X AF X A\ ; (8.20) 



48 



which leads to the result for u. By (8.10) and (8.20) 



( J^) 2 = [4A 2 FA + 4AFA 2 - 8AFAFA\ . (8.22) 

For comparison we have J^- = 2|_A 3 J; hence we obtain (8.21). 

Moreover by Lemma 8.4, |_AJ contains u(x,t 3 ) = 2Ce~ xA FAFe~ xA B and all its deriva- 
tives. Observe that 

-2— v(x, t 3 ) = -4LAJ = u(x, t 3 ) (8.23) 

belongs to |_AJ and satisfies the identity (8.18); moreover, all the partial derivatives of u also 
belong to the differential ring [ AJ . 

□ 

We now point out some particular solutions which are realised via Lemma 8.3, some of 
which were also noted by Poppe [49]. Let Xj be distinct complex numbers for j = 1, . . . , m, 
such that ?R,Xj > 0, and let H = span{x J e~ A£a: : j = 0, . . . , rii — 1; i = 1, . . . , m}, which we 
view as a subspace of L 2 (0, oo), and let A = — ^ on H . 

Corollary 8.4 (Solitons) (i) Then (e~ sA ) se -R defines a Co group of operators on H such that 
\\e- sA \\ < 1 for s > 0, and (f>(x;t 3 ) satisfies §^ = 20, and u(x;t 3 ) G C(x, t 3 , e^ x , e~ 2A ?* 3 ) 
satisfies KdV. 

(ii) In particular, suppose that A has distinct and simple eigenvalues, and that B = 
(bj) G C nxl and C = (Cj) G C lxn . Then 

N - \ - \ 

det(J + nR x ) = » m E II bjCje-***-**" I] ( 8 - 24 ) 

m=0 CT C{l,...,JV},HCT=m jeer j,kea:j^k j k 

Proof, (i) The group e~ sA operates as translations e~ sA f(x) = f(x + s), and hence e~ sA 
is a strict contraction on the finite dimensional space H = C n for s > 0. In effect, we have 
returned to the setting of Proposition 4.4. The generator is —A = d/dx, and can introduce 
A 3 = —d 3 /dx 3 and the group e~ taA which is associated with the linearized Korteweg-de 
Vries equation. By Lemma 8.3, u satisfies the KdV equation (8.23), and by Theorem 4.2, u is 
rational in the basic variables. 

(ii) Apply Proposition 4.4 and Lemma 8.3. 

□ 

Let H = L 2 (— oo, oo) and as in section 6, we can take Af{x) = —f'(x) and we note that e~ tA 
is the Airy group 

i r°° 



e- tA f(y) = ^J /(fle - ** + ^ df. (8.25) 
Then with g G ©(A 4 ) we choose Bo : a H- g(y)ct and Cq : f (->■ /(0), and let 
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Corollary 8.5 (Non solitons) Let g G V(A 4 ) have Y^=oO- + n )l7«l < 1- Then (f>(x; t) satisfies 
(8.16) and u(x;t) satisfies (8.17). 

Proof. By Plancherel's formula we identify V(A 4 ) = {g G L 2 : /^(l + £ 8 )|<K£)| 2d £ < 
so the maps are well defined. By a simple calculation, have 

/>oo 

Rxf(y)= g(y + s)f(s)ds, (/GL 2 (0,cx))) (8.27) 

./ a: 

so in particular Rq is the Hankel integral operator with kernel g(y + s). Hence Ro is unitarily 
equivalent to [7j+fc]°° fc=0 on £ 2 , which by the hypotheses is a trace class operator; likewise, R x 
is trace class. Furthermore, I + R x is invertible, and the inverse F is given by a Neumann 
series. Given these facts, we can apply Lemma 8.3. 

□ 

For m > 4, We can choose g(x) = l^^(x)x m e' x in Corollary 8.4. Whereas the choice 
of g(y) = So is technically inadmissible, the resulting expression 4>(x;t) = t _1 / 3 Ai(— x/(6t) 1 / 3 ) 
does give a solution of (8.16). 

Proposition 8.5 Suppose that CqA 6 : H — >■ C and A 5 Bq : C — >■ H are bounded. 

(i) Then the scattering function </>(x;t 2 ) = Coe -2 * 2 " 4 ~ xA Bq satisfies 

(ii) Let v(x) = T(x, x), so that 

v(x,t) = -C e- xA - t * AB (I + R)- 1 e- xA - t * A ' i Bo. (8.29) 

Then u(x,ts) = |^ satisfies the KdV(5) equation 

du d 5 u d 3 u dud 2 u ,<9« . 

OT5 <7X D ox ox ox z ox 

Proof. We shall prove that 



dv d 5 v d 3 v dv fd 2 v\ 2 fdv\ 6 



16— = ^ + 10 



(O V\ A /ov\° , 



dt 2 dx b dx 3 dx 
The basic identities required follow from (8.20), namely 

= 16LA 4 J - QA[A 3 F X A + AF X A 3 \ - 96[A 2 F X A 2 \ 

+ 112[A 2 F X AF X A + AF X A 2 F X A + AF X AF X A 2 \ - 384[AF X AF X AF X A\ ; (8.32) 

|^ = 32 [A 5 J - 160 [A 4 F X A + AF X A A \ - 320[A 3 F X A 2 + A 2 F X A 3 \ 

+ 640 [A 3 F X AF X A + AF X A 3 F X A + AF X AF X A 3 \ 
+ 960[A 2 F X A 2 F X A + A 2 F X AF X A 2 + AF X A 2 F X A 2 \ 

- 1920[A 2 F X AF X AF X A + AF X A 2 F X AF X A + AF X AF X A 2 F X A + AF X AF X AF X A 2 \ 
+ 3840[AF X AF X AF X AF X A\. (8.33) 
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Using these, one checks that (8.30) holds. 

□ 

9. The Kadomtsev— Petviashvili equations 

In Proposition 5.4, we proved that r functions associated with Schlesinger's equations 
extend to the Jacobians of hyperelliptic curves. In this section, we give a criterion for a r 
function from a linear system to be the restriction of a theta function on the Jacobian of a 
complete algebraic curve. Schottky [45, 54] asked how one can characterize the 9 functions that 
arise from compact Riemann surfaces amongst all the possible functions 6{x | O) on Abelian 
varieties as in (4.1). Shiota proved that the Kadomtsev-Petviashvili system of differential 
equations [54, 62] characterize those r functions that arise from complete algebraic curves. 

Definition (Theta functions) For < g < oo, let A be a lattice in C 9 such that C 9 /A 
is a complex torus, which is compact for the quotient topology. A quotient of nonzero 
entire functions on C 9 is said to be a theta function if there exists a family of linear maps 
C 9 -> C : z L 7 (z) for 7 G A and a map J : A — > C such that 6(^+7) = e 2 ™(M*)+^(7))e(z) 
for all 7 G A and z G C 9 . If Q is a quadratic form on C g , if) : C 9 — >■ C is a linear functional and 
c G C", then e 27rl (<2(^)+V>0)+ c ) gives a trivial theta function, namely a Gaussian. Evidently 
the product of theta functions is again a theta function. In this section, we consider only entire 
theta functions. 

Definition (Riemann's theta function) Suppose that f2o and ^1 are real symmetric j x j 
matrices with Qi positive definite, and let Q, = Qo + ifii; then let A = Z 9 + QJL 9 be a lattice 
in C g . Then 

e( x \n)= e 2vrimtx+7iimtnm (9.1) 

is Riemann's theta function for the Abelian variety X = C 9 /A. Let co G C have Qui > 0; then 
Jacobi's elliptic theta function 9± for the torus C/(Z + 00Z) is 

00 

1 (x\w)=i ^ (_ 1 )n e (2n-l) 7 r^+(n+l/2) 2 7 r^_ ^ 
n= — oo 

Remark Gaussian functions give rise to theta functions which are regarded as trivial. We 
can realise these from linear systems as follows. Given any N < 00 and a positive definite real 
symmetric matrix Q with inverse Q~ l , we introduce a linear system with state space L 2 (R. N ), 
with state variables (x, t) = (x, t\, . . . , fjv-i) and £ = (£0, • ■ • , £jv-i), by 

B:C^H: a{2 N n N det Q)~ 1/4 exp^Q' 1 ^, £)/ 4 ) 

JV-l 

U(t)e- xA U(t) :H^H: /(£) ^ exp(-ix£ - i ^ 

3=1 

C.H^C: f„ / R ,/(O«p(- -(»)/4) (2 ^^-; /< . 
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For consistency with the theory of this section, we define this the tau function of the Gaussian 
linear system to be 



t (x, t) = CU(t)e~ xA U(t)B = exp(-Q((x, t), (x, i))/2). 

For any compact Riemann surface £ of genus g, one can define a homology basis and a 
(/-dimensional space of Abelian differentials of the first kind. Then one defines a corresponding 
lattice A of periods and a Jacobi variety J = C 9 /A with a period matrix f2, and hence 
Riemann's theta function 8(x | Q) by (9.1). When g = 1, any Abelian variety X = C/(Z2wi + 
Z2cji) is the Jacobi variety of some elliptic curve, and in section 11 we realise 9\ as the theta 
function of a periodic linear systems. 

More generally, suppose that r(t) = det(/ + e z ^J=° 3 E) with t G T°° is the tau func- 
tion that arises from a periodic linear system (—A,U(t)B,CU(t);U(t)EU(t)) as in Theorem 
7.2. Given a linear map p : C 9 — > C°° of rank g such that p(ej) G Z°° has only finitely many 
non-zero entries with respect to the standard bases, then p t : C°° — > C 9 satisfies p*(Z°°) C Z 9 . 
Then rop: C 9 — >■ C is entire and periodic with respect to Z 9 . 

Proposition 9.2 (Shiota) Suppose further that r o p(t) = 9(t \ Q), where 6 {t \ ft) is Riemann's 
theta function for some Abelian variety X = C 9 /A of dimension g; let Q(x, y, s) be a quadratic 
formJet k, 7, 8, ( G C 9 with n / 0, and for 

a(x, y, s- C) = e Q{x ' y ' s) 6(Kx + 7 y + 6s + ( \ Q), (9.3) 

let u(x, y, s; () = —2^- loga(x, y, s; (). Then the following two conditions are equivalent: 

(i) the 9 divisor is irreducible, and u satisfies the KP equation 

ox\ox 6 ox ox os J oy z 

for all C G C 9 , in the special case of a = — 1, A = 0, /3 = 1; 

(ii) X is isomorphic to the Jacobian variety of a complete algebraic curve. 

Proof. See the papers of Shiota [54] and Mulase [45] . 

□ 

Proposition 9.2 concerns entire r functions. In section 11, we consider the more topic of 
functions that arise from quotients of Abelian 6 functions, especially the 6 quotients on the 
Jacobian of a compact Riemann surface. 

In the case of (2, 2)-admissible linear systems, one can produce solutions to the KP 
equations via scattering functions and the Gelfand-Levitan equation. Zakharov and Shabat 
[59] considered the associated scattering function which satisfies the linear KP equations 



Oi-p, h 



dt dx 3 dz 3 \dx dz ) 

and 



M-5T + -57)=° ( 9 - 5 ) 



dy dx 2 dz 2 
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Definition The appropriate version of the Gelfand-Levitan equation for the linear KP equa- 
tions is 

poo 

y(x,z;y;t) + K(x,z;y;t)+ / K(x, s; y; t)*(s, z; y; t)ds = {x < z). (9.7) 

J X 

For a solution K(x, z; y; t), define the potential by u(x; y; t) = —2-^K(x, x; y; t), so that u ^ 
is the scattering transform. 

Theorem 9.3 (i) Let (—Ai,B,C) and (—A2,B,C) be (2,2) admissible linear systems as in 
Lemma 2.2 with input and output spaces C, where A\ and A^ are bounded on a common state 
space H. Let C(y;t) = Ce^+^^-y^^ and B(y;t) = e t(A 3 2 +\A 2 )/a+ y Al/p B _ 

(i) Then ^(x,z;y;t) = C(y;t)e~ xAl e~ zA2 B(y;t) satisfies the scattering equations (9.7) 
and (9.8) for the KP equation. 

(ii) With S x = e- A * s B(y;t)C(y;t)e~ A i s ds, let 



K(x, z- y; t) = -C(y; t)e~ x ^ (I + S^e'^ B(y- 1). 
Then K(x, z; y; t) satisfies the integral equation (9.7) and 

K(x, x; y;t) = ^- log det(J + S x ). 
dx 

Moreover, there exists xq such that K(x, z; y; t) satisfies 

d 2 K d 2 K BK . , 

~dx 2 ~ ~ ~dz 2 ~ + ~dy = X ' V] l ' ^' Z] m ^ i x o < x < z ) 



(9.8) 



(9.9) 



(9.10) 



(Hi) Suppose furthermore that A\ = A2 = A. Then the scattering function satisfies 
ty(x,z;y;t) = <p(x + z;t) with 4>(x;t) = Ce 2t( - A + XA )/ a - xA B, and the potential is given in 
terms of a Hankel determinant by 



u(x;t) = -2— logdet(/ + r 0(a!) ). 



(9.11) 



Proof, (i) Since the operators are all bounded, the functions are differentiable and one can 
verify the differential equations, without assuming that A\ and A 2 commute, 
(ii) The linear system 





'-At 




B 




-C 


( 


-^2. 




.0 B(y;t)_ 




C(y;t) 



is (2, 2) admissible and by Lemma 2.2 the resolvent operator 



J~e- sA >B(y;t)C(y;t)e- sA ids 



r 

Jx 



-sAi 



BCe 



-sAn 



ds 



(9.12) 



(9.13) 



53 



is trace class as in Proposition 2.5. Hence one can verify the integral equation as in the proof of 
Proposition 2.4. One then verifies the determinant identity (9.11), which involves the bottom 
left entry S x satisfying the asymmetric Lyapunov equation 

^S x = -A 2 S X - S X A X = -e- xA *B(y; t)C(y; t)e~ xA ^ (x > 0). (9.14) 

The solution of the integral equation is unique for large enough x since ||e _:Eyil || — > 
and He" 2 ^ 2 )! — > exponentially fast as x —> oo; hence \I/(x; z; y; t) — > exponentially fast as 
x — > oo. Using the scattering equation (9.8), one shows by differentiating (9.9) repeatedly that 
-§^K (x, z; y; t) - -§^K(x, z; y; t) - f3-^K(x, z; y; t) and u(x; y; t)K(x, z; y; t) both satisfy (9.9) 
multiplied by u(x; y; t), and so by uniqueness are equal. 

(iii) The determinant det(I + S x ) of (ii) is not generally a tau function in the strict 
sense of Theorem 1.1, but in the following case where A\ = A 2 we do obtain a genuine 
tau function from a linear system with input and output space C. Here the linear system 
{—A,e t<yA + XA y° i B,Ce t ( A ' +XA ^ a ) has scattering function <j)(x;t), and the resolvent operator 
simplifies to 

/>oo 

R x = e^ A3+XA ^ a e- sA BCe' sA e t( - A3+XA ^ a ds, (9.15) 

J X 

hence det(J + S x ) = det(I + T^^) by Proposition 2.4. 

□ 

A solution to KP that is given by a rational curve with only ordinary double points is 
known as a soliton solution; see [45] . By introducing a suitable linear system in discrete time, 
we can recover several examples of soliton solutions u(x;y;t) of KP that are presented by 
Zakharov and Shabat [62]. In particular, suppose that H has dimension n and A G B(iJ); we 
do not need any further assumptions on A. Let B be the controllability matrix 

B=[B AB A 2 B ... A n ~ 1 B] (9.16) 

and with C(y;t) = Ce yA2 ~ tA3 , let C(y;t) be the observability matrix 

C(y;t) = column [C(y;t) C(y;t)A ... C(j/; tJA"" 1 ] ; (9.17) 

Then we can take A 2 = and Ai = A in the context of Theorem 8.3, and simplify the 
expression in () by taking S(x;y;t) = BC(y;t)e~ xA , which satisfies a one-sided version of 
(4.36), namely |f = -SA. 

Corollary 9.4 Suppose that H has dimension n. Then det(J + /j,S(x; y; t)) is a polynomial in 
/i of degree less than or equal to n, and the coefficient of jjL n is a Hankel determinant r n (x; y; t), 
where u(x; y; t) = —2^- logr n (x; y; t) satisfies KP. 

Proof. Let L = [C{y; t)A^ k - 2 e- xA B\ n j k=1 and observe that 

det(/ + iiS{x; y; t)) = det(/ + n.BC(y; t)e~ xA ) 
= det(I + fiC(y;t)e- xA B) 

= det(/ + /iL). (9.18) 
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Since T is an n x n matrix, this gives a polynomial of degree less than or equal to n, with 
equality only if the ranks of B and C(y;t) are both equal to n, in which case the leading 
coefficient is 

r n (x;y;t) = det[C(y;t)Ai+ k ~ 2 B}l k=1 . (9.19) 

Alternatively, we can view r n as the Wronskian with respect to x of the functions ^k{x; y; t) = 
C(y;t)A k ~ 1 e~ xA B,. so by a result of Freeman and Nimmo [24], r n satisfies Hirota's bilinear 
form of the KP equation 



dr n \- 
dy) 



and hence u satisfies KP. 



d 2 



-4tv 



dx 2 
( d 2 r n \ ■ 
V dxdt J 



+ 3r n 



dy 2 
dx 



dt 



dr n \ ( dr n 
dx 



d 4 r n 
dx 4 



= 



(9.20) 



□ 

The KP differential equations reduce to KdV equations in specific cases, and the KdV 
hierarchy is specifically associated with hyperelliptic curves. By considering the specific form 
of Hankel operators, we deduce that hyperelliptic curves give the theta functions that are most 
naturally associated with Hankel determinants as in Proposition 2.4. 

Following Krichever and Novikov [38], we consider the operators 



Lx = — 

dx 



L 3 = 



d_ 

di 





u — k 



k 2 



I - 9 
dy 



-k 







-A: 



1 du 
4 dx 

ku y?_ _i_ 1 d 2 

2 9 ~r 4 



-k 2 

_ 1 du 
4 dx 



note that k \— > Lj is a polynomial for j 



2 1 4 dx 

1, 2, 3, and that trace(Lj 



for j 



that (f)(x;t) = Ce XAt+2A ^"e^B, and let u(x;t) satisfy 

1 d 3 u 



du 
~dt 



4dx 3 



3 du 
2 U dx~- 



(9.21) 

(9.22) 
1, 2. Suppose 

(9.23) 



Then ^(x, z; y; t) = (f)(x + z;t) gives a solution of the scattering equations for KP which does 
not depend upon y. Likewise, u gives a solution of KP which does not depend upon y, and 
the operators L\,Li and L3 commute. 

The KP equation (9.4) is the first of a hierarchy of differential equations. In the matrix 
model for two dimensional quantum gravity, a r function of the KP hierarchy appears as the 
square root of a partition function. 



10. The Baker— Akhiezer function 

Krichever [35] defines the Baker- Akhiezer function ip{x] p) for spectral data (£, po, D, xq) 
as follows. Given a compact Riemann surface £ of genus g with distinguished point po = 00, 
let the local parameter at po be z = A -1 and suppose that p (->• ^(x;p) is meromorphic 
on £ \ {po} with poles 71, . . . , j g which are independent of x = (x±, . . . , xjv); suppose also 
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that there exist £j(x) for j = 1,2, . . ., which are holomorphic for x in a neighbourhood of 
(xi t0 , . . . , Xjv,o)j suc h that there is a convergent series expansion 

Mx, p) = (l + £ ^r) exp(j> fc - x kfi )X k ) (p p ) (10.1) 

i=i fc=i 

see [35, 54]. Given x, xo and the positive divisor D = Y^j=i ^j^u > there exists a unique ^(x, p) 
on £\{p } with these properties, up to a scalar multiple, and one can construct such a function 
from quotients of Riemann's theta functions on the Jacobi variety for £. Krichever showed that 
potentials which are finite gap or algebro-geometric on a complete algebraic curve give rise to 
an algebraic family in the sense of Lemma 5.3, and he constructed the meromorphic matrix 
function W from the Baker-Akhiezer function and its derivatives. To deal with commuting 
families of differential operators of rank greater than one, he introduces matricial ip(x, A) and 
matricial divisors in [36]. 

In contrast, we define our ipBA for linear systems, irrespective of whether there exists 
a suitable £, and then seek to recover the properties of Krichever 's definition. We start by 
looking at admissible linear systems, and then consider periodic linear systems. Then we 
introduce more parameters into the discussion, and consider iPba(x, t, A) with t in an infinite 
dimensional torus. Finally, we recover the case in which the linear system is associated with a 
compact algebraic spectral curve. 

In this section we consider scattering functions and the multiplication rule for tau functions 
which is analogous to the addition rule for positive divisors divisors on an algebraic curve. The 
multiplication B (->■ (XI — A)(XI + A)~ x B is associated with adding the divisor associated with 
a pole on the spectral curve. There is a consequent formula for addition of divisors, which 
Ercolani and McKean [21] credit to Darboux, as in Proposition 4.2. To subtract divisors, one 
considers the quotient of tau functions, as in the definition of the Baker-Akhiezer function. 

Definition (Baker-Akhiezer function) Given an admissible linear system Sqo = (—A, B, C) 
with tau function T ao (x) = det(7 + r^ (x) ) as in Proposition 2.2, we introduce 

E A = (-A, (XI + A)(XI- A)- l B,C^j (3?A>0) (10.2) 
with tau function t\(x), and the Baker-Akhiezer function 

i> BA (x; A) = exp(Ax) (10.3) 

Too [Xj 

Let C^°(R; R) denote the space of infinitely differentiable functions / : R — >■ R such that 
\x\ j f^(x) — >■ as x — >■ ±oo, and suppose that u G C^°(R; R). Then with A = k 2 , let s(k) be 
the scattering matrix, which depends analytically upon k, and let S2i(k) be the bottom left 
entry, which satisfies s 2 i G C^°(R;R) and S2i(k) = S2i(—k), so that 

1 f°° 

<t>{x) = — e lkx s 2l (k)dk (10.4) 

IK J-oo 
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gives a real function. Dyson inverted the scattering map q (->■ s 2 i by the formula (1.12). 

10.1 Admissible linear systems in scattering As in [9, Theorem 4.2] and [21, p. 486] 
we can introduce a linear system and Hankel determinant to realise scattering functions. The 
following formulas are similar, but slightly different from those in [21]. Let H = L 2 (R;C) 
and let b\, b 2 : R — >■ C be smooth functions of compact support such that bi(—k) = bi(k), 
b 2 (-k) = b 2 (k) and \h(k)\ = \b 2 (k)\ for all k G R, and let 

B:C^H:a^ h(k)a; 

e~ xA : H -> H : f(k) ^ e ixk f(k); 

1 f°° 

C: H -)> C :f(k) >-> — / f(k)b 2 {k)dk. (10.5) 
2vr J_ 00 

The potential u is in C^°(R;R), and we assume that there are no bound states, so we are 
in the scattering case of Schrodinger's equation. Then (—A, B, C) has scattering function 
<j)(x) = e txk b(k)dk/2n, while Sj K = (—A, (ihiI—A)(iK,I+A)~ 1 B, C) has scattering function 
4>in(x) = e lxk b(k)(K+k)(K—k)~ 1 dk/27r, which is unambiguously defined for real k since the 

Hilbert transform is bounded on H; the corresponding potential is Ui K (x) = — 2^- logTj K (x). 

The Bloch spectrum of Schrodinger's operator consists of those A G C such that there 
exists a nontrivial and bounded solution / of —f" + uf = Xf. In this case, the Bloch spectrum 
is a double cover of [0, 00) given by ±k \— > k 2 , where ±fc is associated with the unique fi K (x, ±k) 
such that —f" K (x, ±k) + Ui K (x)fi K (x, ±k) = k 2 fi R (x, ±k) and fi K (x, ±k) — e ±lkx — >■ as x — >■ 
±00. The point k is associated with the function (k + K)/(k — k) which has a simple pole at n. 

Proposition 10.1 (i) Suppose that the operator G : L 2 (0, 00) — > L 2 (0, 00) defined by Gf(x) = 
f(x) + T(x, y)f(y)dy is invertible. Then there is a gauge transformation 

G-\-d 2 /dx 2 + u)G = -d 2 /dx 2 . (10.6) 

(ii) The multiplication rule 

s 2 i(k)^^s 21 (k) (10.7) 
k — k 

is equivalent to the addition rule u(x) t-> Ui K (x) for potentials as in 

d 2 

-2— log ^ba(x, in) = Uoo(x) - u iK (x). (10.8) 

(Hi) The Baker-Akhiezer function is given as a series of Fredholm determinants and sat- 
isfies i^ba(x, ik) — e %kx — >■ as x — >■ 00 and 

—ifj' BA (x,ik) + u(x)'4!BA(x,ik) = k 2 i^BA(x,ik) (x G R). (10.9) 

Proof, (i) The operators —d 2 /dx 2 and —d 2 /dx 2 +u are essentially self-adjoint on C^°(0, 00), so 
the identity foo(x, k) = G(e lxk ) for the eigenfunctions shows that G gives a similarity between 
operators on L 2 (0, 00). 
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(ii) We can express the difference of the potentials for the systems as 

u 00 (x)-u iK (x) = -2—\og^-^-, 10.10 

dx 2 T iK (x) 

and then simplify the expressions. 

(iii) With Ti K and the corresponding potential Ui K (x) = — 2^- logTj K (x) defined for the 
linear system T, iK , we introduce 

/>oo 

f tK (x,k) = e lkx + / T iK (x,y)e ik y dy. (10.11) 

J X 

By repeated integration by parts, one verifies —f-' K (x,±k) + Ui K (x)fi K (x,±k) = k 2 fi K (x,±k) 
and fi K (x, ±k) — e ±lkx — >■ as x — >■ ±oo. In particular, with in = oo we can express 

/>oo 

f^x, k) = e ikx - Ce- xA (I + R x y l / e - yA Be iky dy 

J X 

= e ikx ((l + Ce~ xA {I + R^-^ikl - A)~ 1 e~ xA B) 

= e ikx det(J + (iJfeJ - AJ-^-^SCe-^C/ + i^) -1 ) (10.12) 



where we have used a simple identity for rank-one operators, hence 

(ikI-A)-\ 
det(I + R x ) 



/-(x, t) = e<fa jg!(£+^±(!^i^gg^), (10 . 13 ) 



and we can finish by using Lyapunov's equation 

f (xk)- c ik X det(I + R x -(ikI-A)- 1 K) 

/oolX,ftj " e det(I + R x ) ' 1 j 

where the determinant on the numerator is 

det(J + R x + (zfcJ - A) -1 (AR X + R X A)) = det(l + R x (ikl + - A) -1 ). (10.15) 

As in classical Fredholm theory [55, p. 70], we can expand the determinant on the denom- 
inator of (10.14). As in [10], the operator R x on L 2 (0,oo) is represented by the continuous 
kernel 

R * M= -2ni(s + t) ' (10 ' 16) 
so we have a Cauchy determinant 

-e^6i(sj)&2(s^e ia:s< 



R x [ Sl ■■ " n ) = det 



Si ... s 



(10.17) 

3,1=1 



-2ni(sj + s^) 
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which is rapidly decreasing as s, t — >■ oo. Hence we have a convergent series 

r 00 (x) = det(I + R x ) = f]- 1 f ... f R x ( Sl Sn )d Sl ...ds r 



n=0 n - J J (0,oo) 

Likewise, we can represent 



(10.18) 



so we can expand det(J + R x (ikl + A){ikl — A)~ l ) as in (10.16). To deal with convergence, 
we can take k = ko + in and let r\ — >■ 0. Hence the Baker- Akhiezer function foo(x; k) may be 
expressed as in (10.1). 

□ 

The pole divisor z\(t) is determined by {z n (t) : ipBA(zn(t), t; A) = 0} and is associated 
with the potential u\(x;t) = — 2-^ log T\(x,t). In this section, we introduce dynamical sys- 
tems on 7^° such that u\(x,t) undergoes the nonlinear evolution associated with the KdV 
hierarchy. To obtain KdV(2n + 1), we vary t2n+i while fixing £2.7+1 for j ^ n. 

Definition (Torus) Let (tqo) = {p n : n = 1, 2, . . .} be the zeros of Too and O n be the real oval 
in U^ £ (_ 00 00 ](ta) that is based upon p n . Then let 7^° = Il^Li On an d consider z\ = {z n : 
n = 1, 2, . . .} = (t\) with z n E O n . Then z\ G 7r° is the pole divisor of i^ba(x, A) in the 
infinite real torus 7^°. 

Proposition 10.2 (i) The Baker-Akhiezer function A t-> 4>ba(x,X) is holomorphic on C \ 
Spec(A). If E has finite rank, then A t-> i/j B a(x, A) is meromorphic on C with the only possible 
poles being on the spectrum of A. 

(ii) i/jba(x; A) belongs to a Liouvillian extension of the field of fractions of Aq and satisfies, 
in the notation of Proposition 2.2, 

i) BA (x,\) = e Xx det(l- J T(x,y)e x(y - x) dy^ (9&A < 0). (10.20) 

(ii) Suppose that S is a block diagonal direct sum (BfLi^j, where £j is a periodic linear 
system with Tj as in Theorem 7.1. Then 

i> BA (x,X) = e Xx JJdetfl- / T j (x,y)e My - x) dy\ (»A < 0). (10.21) 

j=l J x 

(Hi) Suppose that B and C have rank one. Then 

-ij' BA (x, A) + u(x)i/>ba(x, A) = -\ 2 ^ba{x, A). (10.22) 

(iv) If t\ has only simple zeros, then each zero of i^ba(z,X) in (t\) processes in a real 
oval based at a pole of iPba(z, A) in (tqo) as A describes (—00, 00). The pole divisor defines a 
map Sa ' ^ z\ from the periodic linear system to the real torus 7^°. 
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Proof, (i) The first part follows from Theorem 7.2. Suppose then that E has finite rank, and 
note that (XI + A) (XI — A)~ l E is a rational function with values in the space of operators 
on a finite-dimensional Hilbert space. Hence the determinant t\ is meromorphic as a function 
of A on P , and the only possible poles are at the A G Spec^). Hence A i— >■ iPba(x,X) is 
meromorphic on C and has poles independent of x, as in Krichever's definition of the Baker- 
Akhiezer function. 

(ii) We have u(x, A) = — 2(logTA)" in A by Theorem 7.2, hence ip'^ A (x, ^) belongs to A ; 
we integrate this to obtain 4>ba in some Liouville extension. By some simple manipulations, 
we have 

det (7 + R X (XI + A)(XI - A)- 1 ^ = det(I + R x ) det (i + (XI - A)-\AR X + R X A)(I + 

(10.23) 

where AR X + R X A = e~ xA BCe~ xA , and hence 



det(l + R x (XI + A)(XI-A)- 1 ) 

'- = det(/ + Ce--(/ + fl.)-(A, - A)~^b) , (10.24) 

and e x ( y ~ x ^e~ yA = —(XI — A)~ 1 e~ xA , which leads to the stated identity. Moreover, the 
right-hand side is analytic in A when |A| > \\A\\, and i^ba(x,X) = e Xx (l + 0(A~ 1 )) as |A| — > oo. 
By the proof of Theorem 7.2, 4-^ logics a(x; A) belongs to A . 

(iii) This follows immediately from (i). 

(iv) We reduce to the case of the admissible linear system (—A — eI,B,C), which has 
input and output space C, as in Proposition 2.5. For e > 0, let R x £ ^ = e~ 2ex e~ xA Ee~ xA , so 
that R x 6 ^ -> exponentially fast as x — >■ oo, and R x 6 ^ satisfies the Lyapunov equations 

-^R^ = (A + eI)B& + R X £ \A + el), (10.25) 

with 

-^-R x £) \ x=0 = BC + 2eE. (10.26) 

ax 

Since BC and E are trace class, we can introduce t&(x) = det (I + R x 6 ^) and 

T ( x e) (x) = det(I + R { x £) (XI + el + A)(XI - el - A)- 1 ), (10.27) 

whenever A— £ is in the resolvent set of A; likewise we can introduce (x) = —2^- log t^} (x). 
Now the Baker- Akhiezer function 

f (e) (x,k) = e ikx ^^- (10.28) 

satisfies 

- ^f {£) (x) + « (£) (x)f^ (x, k) = k 2 f^ (x, A;); (10.29) 
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letting e — >■ 0, we obtain 



as required. 

(v) Clearly the poles of i^ba(z, A) occur at the zeros of Too(z), and hence form the set 
(too), for all A. The zeros of i^ba(z, A) form the set (t\), which does vary with A. The subset 
{(XI — A) (XI + A)~ l : A G R} of B(H) is compact in the norm topology since the spectrum 
of A is separated from R; hence t\ gives a compact family of holomorphic functions for the 
topology of uniform convergence on compact sets, with T- 00 (z) = Too(z). For each bounded 
open subset f2 of C, the set {zGfi: t~\(z) = 0} has a uniformly bounded number of terms for 
— oo < A < oo, by Jensen's formula and the Proposition 7.4. Each zero depends continuously 
upon A by the inverse function theorem, and describes an oval for — oo < A < oo. 

□ 

Given a tau function from a periodic linear system (—A, B, C; E) where the input and output 
space are both C, we consider the conditions under which r arises from the theta functions 
on a compact algebraic curve. First we consider families of linear systems as in Theorem 10.2, 
with common A, which are parametrized by A G P 1 \Spec(v4) and time parameters (t3,t$, . . .), 
giving tau functions T\(x,t). Initially x and £2.7+1 are real, and T\(x,t) is n periodic in each 
variable, hence r\(x,t) gives a periodic function on the infinite real torus H°° / ttZ 00 ; then 
we extend to complex x and £2.7+1, so that t\(x, t) is entire. By forming quotients of such 
functions, we aim to realise typical tau functions. 
To introduce the required linear systems, we let 

T = {(x,t 3 ,t B , . . .) G R°° : lim sup \t 2j+ i\ 1/j = 0} (10.28) 
which gives an Abelian group under addition, and for (x, t) G T, let 

U(t) = exp(— Y^jLi hj+iA 2 ^ 1 ), which gives a multi parameter group of operators such that 
U(s + t) = U(s)U(t). Then we replace £00 (0) = {-A, B, C; E) of Theorem 10.2 by 

S A (t) = {-A, (XI + A)(XI - A)- l U(t)B, CU(t), (XI + A)(XI - Ay l U(t)EU(t)^j (10.29) 

for A G P 1 \ Spec(^4). Each £x(£) gives a space Ao(t, A) of potentials as in Theorem 10.2(iii), 
while A is a spectral parameter as in Proposition 10.5. Let (A ,d/dx) be the differential ring 
generated by S as in Theorem 10.3(h), and let (A 00 ,d/dx,d/dt 2 j+i) be the differential ring 
generated by all the T,\(t); then A C Aoo, and the inclusion splits by mapping £2^+1 ^ for 
all j = 1,2,.... 

Definition (Baker-Akhiezer function) We define the quotient 

xX + h j+ i\ 2j+1 ) —± , L (10.30) 

' Tco \X, t3, *5 , . . .) 
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to be the Baker-Akhiezer function of the periodic linear system (—A, B, C; E) under U(t). 

This definition is consistent with earlier in this section, and with Shiota's definition [54]. 
However, we cannot expect a precise analogue of Proposition 10.2, which expresses eigenfunc- 
tions in terms of i^ba- 

Lemma 10.3 (i) The scattering function §\(x,y) = CU(t)e- xA (XI + A)(XI - A^UifyB for 
satisfies 

£^* A(:c ' t)+ s^* A( ' T '' )=o - (10 - 31) 

(ii) t\(x, t) is holomorphic for (x, t, A) G C x C°° x (P 1 \ Spec(A)), where C = C/nZ is the 
complex cylinder. 

(Hi) X i y if)BA(x,t, X) is meromorphic on C \ Spec(A), while (x,t) t-> if)BA(x,t, A) is 
meromorphic and quasiperiodic with respect to the lattice tvZ°° in C°°. 

Proof, (i) Since U (t) is actually analytic in each tj this is a straightforward computation, 
(ii) First we observe that 

t x (x, t) = det (7 + (XI + A)(XI - Ay 1 U(2t)e- 2xA E^j . (10.32) 

as in Lemma 8.1. Hence A (->■ r\(x,t) is holomorphic on P 1 \ Spec(A), and (x,t) (->■ r\(x,t) is 
entire in each variable since A is bounded. The spectrum of A 2 ^ 1 is contained in 
{-iN 2 '* 1 , -i(N - l) 2j+1 , . . .,iN 2 i +1 }, so e 27rA - 2j+1 = I, and by Theorem 10.2 r A (x + ir,t) = 
T\(x, t); likewise t\(x, t) is unchanged by adding n to tj] so t\(x, t) is periodic with respect to 
vrZ 00 in C°°. 

(hi) The function YlJLi ^2j+iA 2j+1 is entire by the choice of (x, t) G T, so A (->■ ^ B a(x, t, A) 
is holomorphic on C \ Spec(A). With (ej)°^ the standard unit vector basis in T°°, we 
deduce from (ii) that ^BA^^t + irej,X) = e 27Tx2j+1 r il)BA(x,t,X), and (x,t) i-» ^BAix.t.X) is 
meromorphic. 

□ 

For any periodic linear system as in Theorem, we let 

t] A; x ) = r ffl "Vfc(x,t;A;x ) l N ^ ^ = (1Q 33) 

L ij,k=i 

be the Wronskian matrix associated with some collection of functions, and suppose that $ 
is invertible; then with Uj (x, t; A; Xq) = Q ®® - the operators Lj = Qt ® - — Uj give a 
commuting family of matrix differential operators. 

Corollary 104.4 Suppose further that there exists a compact Riemann surface P with marked 
points {pi, . . . , p m } and field K of meromorphic functions on P\{pi, . . . , p m } and that W(t, p) 
is a M]\f-valued function on T°° x P such that the entries of p W(t, p) belong to K, and 
that [Lj,W] = 0. Then the family Lj is algebraic and has a spectral curve £, with marked 
points {qi, . . . , q n m}, that gives a finite cover of P. 

Proof. This is immediate from Lemma 5.3. 

□ 
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Consider Hill's equation in the form 



= 0, 



d 

dx 



1 

u-X 



(— OO < X < 00 ) 



(10.34) 



where u is continuous, complex- valued and 7r-periodic on R, and let F\(x) be the fundamental 
solution matrix so that LF\(x) = and F\(0) = 1^- Let W be the monodromy operator 
W : ty(x) i y \I/(x + 7r), which commutes with L. Now H\ = : L\I/ = 0} is a two-dimensional 
complex vector space, hence W\h x may be represented by the monodromy matrix -Fa(7t), and 
the discriminant is A(A) = trace F\(n). Suppose that F\(ti) has eigenvalues p and 1/p, and 
that p + 1/p G [—2, 2]; then p\ = 1 so -Fa(^) is uniformly bounded for real x. Hence the Bloch 
spectrum {A G R : A (A) 2 < 4} consists of those points such that Hill's equation has a pair of 
independent bounded solutions. Each oval O n is associated with a gap in the Bloch spectrum 



Definition The multiplier curve is £ = {(A, p) : p 2 — A(A)yU + 1 = 0}, and potentials are said 
to belong to the same spectral equivalence class if their multiplier curves are equal. 

We now consider how the results of section 7 relate to the notions of Liouville integrability 
and finite gap integration. The results of this section are essentially corollaries of some subtle 
results proved elsewhere, and the most interesting relate to elliptic potentials. 

The solutions of (10.1) turn out to be complicated polynomials in u and its derivatives, as 
one can prove by induction. Nevertheless, we can express a solution g m simply in terms of |_AJ . 
The following proposition is a compilation of known results, and included for completeness. 

Proposition 10.5 Let (—A, B, C; E) be a periodic linear system as in Theorem 7.2. 

(i) The complex vector space spanned by the [A 2171 ^ 1 ] for m = 1,2,... is finite dimen- 
sional. 

(ii) Suppose further that S is commutative, as in Corollary 7.4, and that the operators 
in [A\ all have a common eigenvector. Then there is a homomorphism of differential rings 
p : |_AJ — > C(Q,; B(H)) such that g m = p([A 2m+1 \) gives solutions to the the stationary KdV 
equations. 

(Hi) Suppose that u has finite gap, so that g m = for some m. Then there exists a 
hyperelliptic spectral curve £ is hyperelliptic, with a branch point po be a branch point, a 
meromorphic function A on £, and a pair of distinct points Pj,qj G £ for each point ij G 
Spec(A), all independent of x, such that A t-> i/jba(x,X) is holomorphic on £ \ {pj,qj ■ j = 
0;ij G Spec(A)}. 

Proof, (i) Let m be the minimal polynomial of degree iV for the algebraic operator A. Then 
for each entire function /, either f(A) = or there exists a polynomial r of degree less than 
or equal to iV such that f(A) = r(A). Hence the span of the A 2m ~ x for m = 1,2,... is 
finite-dimensional, and hence its image under |_ . J is also finite-dimensional. 

(ii) By Theorem 7.2(h) and Corollary 7.4, |_AJ is a commutative differential ring of opera- 
tors on H over Q. Let ip be a common eigenvector and let p( [P\ ) = A(P) where [P\ ip = X(P)ip, 
so that p defines homomorphism of commutative rings. By Proposition 3.4. the images 
9m = p([A 2m+1 \) satisfies the recurrence relation for the KdV hierarchy. 



[21]. 
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(iii) From the recurrence relation (3.15), we deduce that g n = for all n > m, so u is finite 
gap and C[A, u, u', u", . . .] = C[A, u, u', u" , . . . , u^ m+1 ^] is a differential ring. Any solution of 
the stationary KdV equations is meromorphic on C [52, 6.10]. Let Ao < Ai < . . . < \2 g be the 
simple zeros of 4 — A(A) 2 = 0, and introduce the spectral curve 

29 

£ = {(z,w) : w 2 = Y[(z- Aj)} U {(oo,oo)}, (10.35) 
i=o } 

Now there exists a solution p(x, A) to Drach's equation (8.2) 

p 2 = ~p{x, X)p"(x, A) + - A p\x, A) 2 + (u(x) + X)p(x, A) 2 (10.36) 

such that p(X) is independent of x and A (->■ p(x, A) is a polynomial, which we factor as 
p(x, A) = nf=i(^ ~~ 1j( x ))- Brezhnev [13] gives the solution 

^ ( x)=ex P (gy \^r)' ( 10 - 37 ) 

where the integral is taken along £. Here u and its derivatives are rational functions on £; see 
[34, 50]. For such a potential u, the functions ip± of (10.37) give locally meromorphic solutions 
to Schrodinger's equation. 

Suppose that £ has genus g > 2, and choose po to be one of the 2g + 2 branch points of 
the holomorphic two sheeted cover £ — > P, and then observe that there exists a meromorphic 
function A on £ such has precisely one pole, namely a double pole at po, and hence has degree 
two (When g = 1, we can use A(p) = p(p — Po))- 

The exponential e xX gives an essential singularity in the variable A for p close to po- Also, 
A i — y (XI + A)(XI — A)~ 1 E is a rational function, with trace class values, and the only possible 
poles are on the spectrum of A; hence p i— >■ t\ gives a holomorphic function, except at finitely 
many points of £, which we list as p-,, qj for ij in the spectrum of A. (These can give essential 
singularities when E has infinite rank.) 

□ 

Definition Say that a periodic linear system (—A, B, C; E) is a Picard system if —ip"+uip(x) = 
X 2 if) has a meromorphic general solution ip for all but finitely many A G C. See [28]. 

Suppose that {—A, B, C; E) is a Picard system. Then by elementary Floquet theory, there 
exists a nontrivial solution tp such that tp(x + tt) = pi^(x) for all x. 

Given u G Kc, one can ask whether u is finite gap, and seek to find the spectral curve. 
Gesztesy and Weikard found a conceptually simple characterization of elliptic potentials that 
are finite gap, namely those that are Picard potentials. In the next section, we realise some 
elliptic potentials u that are finite gap in terms of linear systems. 

11. Linear systems with elliptic potentials 

In this section we produce explicit examples of periodic linear systems such that u is finite 
gap, and the corresponding spectral curve £ is of arbitrary genus. 
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Definition (Elliptic functions) Suppose that A = Z2ui + Z2u> 2 with £5(0*2 /wi) > is a lattice, 
and let T = C/ A be the torus, and C = Z/27rZ the cylinder. A meromorphic function on C is 
elliptic (of the first kind) if it is doubly periodic with respect to A; let KL be the differential 
field of elliptic functions. A meromorphic function is elliptic of the second kind if there exist 
multipliers pj G C such that f(z + 2ujj) = pjf(z); so that / is quasi-periodic with respect to 
the lattice; let Kj- be the field of elliptic functions of the second kind. Also let Kj- be the 
set of elliptic functions of the third kind, namely the meromorphic functions on C that satisfy 
f(z + 2ojj) = e ajZ+bj f(z) for j = 1, 2 and some a,j, bj G C. Let Mc be the differential field of 
27r-periodic meromorphic functions; then Kj- C Kj- C Kj- C Mc, where all there spaces are 
closed under multiplication. See [39]. 

First we shall obtain a representation for the coordinate ring Cj- of regular functions on 
elliptic curve 



T = {(X, Z) : Z 2 = 4(X - ei )(X - e 2 )(X - e 3 )} U {(oo, oo)}. 



(11.1) 



Let 9\ be Jacobi's elliptic theta function, Q\(z) be the entire function 0*(z) = 0\(z) and 
let p be Weierstrass's elliptic function with real constants 63 < e 2 < e±. Then (p') 2 = 
4(p — ei)(p — e2){p — e^) so a typical point on T is (X, Z) = (p, p'); moreover KL = C(p)[p']. 

Definition (Realising elliptic theta functions) (i) We refine the basic construction from [11] 
so as to ensure that the various matrices commute. Let H = 
of column vectors and let 



5^L C be expressed as a space 



J = 







I = 



(11.2) 



then for an elliptic nome < q < 1, we introduce the block diagonal matrices on H with 2x2 
blocks, in which each top left corner is exceptional: 





r(i/2)j 













-il 















j 





... 




2q 2 / 


















j 


... 


B = - 





2q 4 / 
















j ... 








2q 8 / 






■/ 








..." 




--u 















j 





... 




q 2 / 








C = 








j 


... 


Eq = — 





q 4 / 


















j ... 








q 8 / 



(11.3) 



Then, with standing for the Hermitian conjugate of A, we introduce 



A = 


~A 



" 

A o_ 


, B = 


~B 



" 

< 


C = 


"Co 



" 

r t 
°o . 


, E = 


Eo 



" 

4. 



(11.4) 
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Given A G C \ {±z}, we introduce a by (XI — J)(XI + J)^ 1 = /cos 2a — Jsin2a; so the 
effect of multiplying B by (XI — A) (XI + A)^ 1 is equivalent to x \— > x + a. 

The spectral multiplicity function v of A satisfies v(i) = v(—i) = oo, which is the antithesis 
of the hypothesis of Lemma 7.1. 

Proposition 11.1. (i) The hypotheses of Theorem 7.2 are satisfied, so e~ xA Ee~ xA defines a 
trace class operator on H , and det(7 + e~ xA Ee~ xA ) is an elliptic function of the third kind 
which satisfies 

oo 

e x {x)dl{x) =det(l + e- xA Ee- xA )\q\ 1 / 2 f](l-q 2n ) 2 , (11.5) 

n=l 

where 6i(x)9\(x) is entire and nonzero on C \ {jir + ifclogq : j,k G Z}. 

(ii) Let S = J£c[I , A, B,C, F]. Then S is a commutative and Noetherian ring of block 
diagonal matrices with entries from Kc; furthermore, S is a complex differential ring for 
(-A,B,C) on C/4ttZ. 

(Hi) The potential u(x) = — 4trace [A\ is the elliptic function 

u(x) = Ap(x) - 4ei - 2(log^i^)"(l/2). (11.6) 

(iv) Then u(x, t) = u(x — ct) gives the general travelling wave solution of the Korteweg-de 
Vries equation 

d 3 u du du . 
dx 3 dx dt 

that has speed c = 4(e x + e 2 + e 3 ) - 3d - (3/2)(log0i^)"(l/2). 

(v) Let A = span c {l, p^\x) : j = 0,1,2,...} and A be as in Lemma 3.2. Then 
A = C [T] , and every element of A with zero constant term is the trace of some element of 
A. 

(vi) The scattering function satisfies 

— 8q 2 

Mx) = \sinx. (x G R) (11.8) 

1 — q z 

Proof, (i) The matrix J satisfies the identities e~ xJ = I cosx — Jsinx and det(J— q 2n e~ 2xJ ) = 
(1 — 2q 2n cos2x + q 4n ). We deduce that e~ xA belongs to S and defines a unitary operator on 
Hilbert space £ 2 ; evidently E is trace class. One can calculate 

oo 

det(/ + e- xA "E e- xA ") = 2i sinx JJ (1 - 2q 2n cos 2x + q 4n ) 

n=l 

101{X) (11.9) 



q 1/4 ir=i(i-q 2n )' 

which reduces to a multiple of Jacobi's function as in [42]. 

Let q = e in ". Then x (x + tt) = -6»i(x) and 6 x (x + 2ttlo) = e- Hx - u ™Q x (x), so Q X Q\ 
is periodic with period tt, and B x (x + 2txu)Q\(x + 2-kuj) = e - 8iix+7TUj) 9 1 (x)ei(x), hence 0i0J 
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is elliptic of the third kind. Using (11.9), one can easily show that the zero set of 6\ is 
{j7T + zfclogq : j, k G Z}, and this coincides with the zero set of 0*. 

(ii) First note that (A 2 + J)(,4 2 + 7/4) = 0, so E = 2~ 1 A~ 1 BC belongs to S. It follows 
directly from Theorem 8.2 that S is a differential ring for (—A,B,C). In this case A is similar 
to — A, so there exists an invertible S such that AS + SA = 0, so the solution to (1.11) is not 
unique. 

Note that the 2x2 matrices satisfy (I + iJe~ xJ )(I — iJe xJ ) = 2isinxl and 

(/ -q 2n e- 2xJ )(I -q 2n e 2xJ ) = (1 - q 2n cos 2x + q 4n )/, (11.10) 

so F is a block diagonal matrix with entries from Kc[7, J]. In terms of t = tanx/2, the n th 
block has determinant 1 + q 4n — 2q 4n (l +i 4 — 6t 2 )/(l + t 2 ) 2 , which has simple zeros and double 
poles for all n. 

(iii) Using the identity (8.1), one checks that 

d 2 

^2 log MI = 2trace[AJ, (H-H) 

then a standard result from elliptic function theory [42, p. 132] gives 

p(x) = -(log0i(x))" + ei + (log0i)"(l/2), (11.12) 

hence the result follows from (11.9). By forming the trace in u = — 4trace[^4J, we are under- 
going a limiting process which takes us from Kc to Me, which includes the elliptic functions. 

(iv) We have the basic differential equation 

p" = 6p 2 - 4(ei + e 2 + e 3 )p + 2(eie 2 + eie 3 + e 2 e 3 ). (11.13) 

One can show that u(x— ct) is a solution by differentiating (11.13) again and then adjusting 
the constants. Conversely, the expression u'" = 3uu' — cu' reduces to 

(u'/4) 2 = 4((u/4f - (c/4)(n/4) 2 + 0(u/4) + 7 ), (11.14) 

where /3 and 7 are constants. By integrating this ordinary differential equation, we obtain 
Weierstrass's function. 

(v) By induction, one can prove that for each n = 0, 1, 2, ... , there exists a polynomial 
q n (X) of degree n + 1 such that p^ 2n ^(x) = q n (X); likewise by induction one can prove that 
there exists a polynomial p n {X) of degree n such that p( 2n+1 )(x) = p n (X)Z. Hence 

span{l, p {j \x) : j = 0, 1, . . . , 2N} C span{X j , X k Z : j = 0, . . . , N + 1; k = 0, . . . , N - 1} 

(11.15) 

and both spaces have dimension 2iV + 2, so we have equality. We deduce that A = {p(X)Z + 
q(X) : p(X),q(X) G C[X]}, which is isomorphic to the ring C[X, Z] modulo the ideal (Z 2 - 
4(X — e\)(X — e2){X — e 3 )); this is an integral domain since Tis irreducible. Hence A = C[7~]. 
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By repeatedly differentiating and using (11.12), we obtain 

p{x) = -trace [A\ + ei + (l/2)(log^^)"(l/2), 

p'(x) = -2trace[A(7 - 2F)A\ (11.16) 

and likewise p^ is the trace of a [Pj\ for some Pj G A for j = 2, 3, . . .. 

(vi) By definition |_-F~ 2 J = Ce~ 2:Ej4 -B. We observe also that C e~ xAo B equals 



-i/cos(x/2) + iJsin(x/2) 

— q 2 Jcosx — q 2 /sinx 

— q 4 Jcosx — q 4 /sinx ... 



(11.17) 



so when we take the trace, we get 

4q 2 

trace C e- xA ° B = -2icos(x/2) -^sinx, (11.18) 

1 — q z 

and we obtain the stated result when we add the complex conjugate to get trace Ce~ xA B. 

□ 

On a compact Riemann surface £, the divisor group D{£) = {5 = Y2j n j^zj : %■ G Z, Zj G 
£} is the free Abelian group generated by the points of £, and the degree of the divisor <5 is 
deg(<5) = n j- We let be the multiplicative group of non zero meromorphic functions on 
£, where we identify / ~ g if / = Xg for some A G C \ {0}. Then by Liouville's theorem, each 
/ G Kje- corresponds to the principal divisor 5(f) = Y2j n j^z.j — X^j m j$Pji where Zj is a zero 
of / of order rij and pj a pole of / of order rrij; moreover, deg(/) = 0. By extension, we can 
consider the notion of a divisor for iPba(z;X) as in Proposition 8.4, with the understanding 
that there are infinitely many zeros and poles on C in a periodic array. In particular, elliptic 
functions of the third kind give rise to divisors on the torus. See [49]. 

We now use the notations r and a to refer to tau functions of periodic linear systems as 
in Theorem 8.2. First consider the group Gc = {t/ct : r, a G Cc} generated by linear systems 
with E of finite rank. Then each rja G K[l may be transformed by the change of variable 
t = tanz to r/d G Kp X and hence gives a divisor 8{r/a) on the Riemann sphere. One can 
check that all divisors of degree zero on the Riemann sphere arise in this way. 

Next consider C/A. The torus T may be identified with the quotient group of divisors of 
degree zero modulo the group of principal divisors, known as the Jacobi variety. We consider 

T M=^' w !Iir'!f7' (n.19) 

which is meromorphic with divisor <5 T ) = J2j — 12k <W on some cen °f the- quotient space 
C/A so deg(r) = n — m. If deg(r) = 0, X^=i( a j ~ bj) G A and a = b = 0, then by 
Abel's theorem r is elliptic of the first kind. In all cases, r is elliptic of the third kind, 
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and u(x) = —2(d 2 /dx 2 )\ogT is elliptic of the first kind with possible poles at the aj, bj and 
congruent points with respect to the lattice A. 

Lemma 11.2 Let ti(x) = e bx+c YTj=i — aj) and To (x) = Ilj=i @i( x ~ bj), and further sup- 
pose that the aj and bj all give distinct points in C/A; next let v(x) = (d/dx) log(ri(x)/ro(x)) 
and w(x) = (d/dx) log(r 1 (x)r (x)). 

(i) Then t(x) = t±(x) / tq(x) is elliptic of the second kind, whereas v(x) and w'(x) + v (x) 2 
are elliptic of the first kind with poles of order less than or equal to one. 

(ii) For each positive divisor (5) on C/A, there exists a periodic linear system with tau 
function tq as in Theorem 7.2 such that (5) equals the divisor of the zeros of To- 

(Hi) Let u be a nonconstant elliptic function on T. Then Ko = C (u) [u'] is a differential 
field, and u = t\/tq where T\ and tq are the r functions of periodic linear systems. 

(iv) Any trivial theta function arises from the quotient of theta functions for Gaussian 
linear systems on R. The effect of multiplying by a trivial theta function r (->■ e~^/ 2 T is to 
take u i— » u + go for some constant go • 

Proof (i) The poles of w and v occur at the zeros of tq(x)t\(x), namely the points aj, bj and 
all the points congruent to these with respect to the lattice A, and the second order poles of 
w' cancel with the second order poles of v(x) 2 . 

(ii) We have obtained a periodic linear system with tau function 9i(x), so by forming 
block diagonal sums, we can obtain a periodic linear system with tau function YYj=i di{x — bj). 

(iii) As in (ii), we can realise u as the u = t\/tq, where r\ and To are r functions of 
periodic linear systems. By a classical theorem, there exists a polynomial P G C[x, y] such 
that P(u,u') = 0, and hence u' is algebraic over C(u), so K is a field, and closed under 
differentiation. In particular, it contains — 2(logTi/Vo) // . 

(iv) See section 9. 



Definition (Differential Galois group [57]) Let U be a fundamental solution matrix of Hill's 
equation 



d 


>" 




o r 






dx 


J. 




u — X 







(11.20) 



with det U(0) = 1 and let PV be the Picard-Vessiot ring over C that is generated by the entries 
of U; then let L be the field of fractions of PV. The differential Galois group DGal(L; C) is 
the set of C-linear field automorphisms of L that commute with d/dx. 

Consider (11.20), where u is elliptic. Then u is the quotient of tau functions from periodic 
linear systems, and may or may not be finite gap. The Treibich-Verdier potentials of the form 



are finite gap if and only if Cj = dj(dj + 1) for some dj G Z for j = 1, . . . , 4, ao G C and the 
the poles satisfy 03 = a\ + 02 and a^ = 0. The corresponding tau function is 



□ 



4 




(11.21) 



4 



r 



(z) = He^z-aj)^^/ 2 eK 3 7 



(11.22) 
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where the exponents are triangular numbers. Whereas one can realise such tau functions from 
periodic linear systems by means of Proposition 11.1, one can likewise produce tau functions 
corresponding to elliptic potentials that are not of the form (11.33). 

We now characterize finite gap elliptic potentials in terms of periodic linear systems. 

Theorem 11.3 Let u be elliptic of the first kind. 

(i) Suppose that (11.20) has a general solution ifj\(x) that is a quotient of r functions 
from periodic linear systems for all but finitely many A G C. Then u is finite gap. 

(ii) Conversely, suppose that u is finite-gap. Then for all but finitely many A G C, (11.20) 
has a solution tp\(x) that is the quotient of tau functions arising from periodic linear systems 
as in Theorem 7.2 and Gaussian linear systems. Also, deg[K^- : K ] is finite. 

Proof, (i) Gesztesy and Weikard [28] considered —ip" + utjj = Xtp for u G K~ ; and showed 
that u is finite gap if and only if u is a Picard potential. If ip is a quotient of r functions, then 
ijj is meromorphic and hence u is a Picard potential. 

(ii) Suppose that (11.20) has a meromorphic solution. Then by a theorem of Picard [28], 
there exists a solution ip that is elliptic of the second kind, hence has the form (11.19) with 
a = and degree zero. By inspecting the differential equation, we see that the only possible 
poles of u are contained in the set {tti, . . . , a n ; b±, . . . , b n }. By Proposition 10.7, each factor 
9±(x — a,j) or 9i(x — bj) arises from the tau function of a periodic linear system, while the 
factor e bx is a quotient of Gaussian tau functions. By [42, p. 96], u' is algebraic over C(it) 
and we have Ko = C(u)[it'] and deg[K^- : Ko] < oo. Let V\ be the solution space of (11.15), 
and observe that DGal(L; Ko) operates on V\ component-wise; in particular, the monodromy 
operators Tj : &(z) ^(z + 2ojj) are commuting operators such that Tj(V\) C V\ for j = 1, 2 
since u is elliptic, so we can take A to be the group generated by T\ and T 2 . Let ^\ ~ \&2 
if <I»i = c* 2 for some constant c G C \ {0}; then let V{ = (V\ \ {0})/ ~. Then with V 
the solution that is elliptic of the second kind, $ = column [ip ip'] G V\, gives a common 
eigenvector T\ty = pi*]/ and T2^> = P2^, so 7^ ~ for all 7 G A; hence $ gives an element of 
(V A * : A). Furthermore, if T\ or T2 has distinct eigenvalues as an operator on V\, then there 
exists a fundamental system of elliptic functions of the second kind, so (V£ : A) is isomorphic 
to P 1 . 

□ 

Remarks 11.4. (i) The fundamental example of a finite gap elliptic differential equation 
is Lame's equation. Let (X,Z) = (p(x),p'(x)) and (Y,W) = (p(y),p'(y)) be points on the 
elliptic curve T = {(X, Z) : Z 2 = 4X 3 - ^ 2 A" - g 3 } U {(00, 00)}, where g\ - TJg\ ^ 0, with 
Klein's invariant J = g\j(g\ — 27g 2 ), and K = C(X)[Z] is the elliptic function field. Then 
Lame's equation is 



dXJ 



1 

£(£+l)X-X 



(11.23) 



and we write ^(X) = ip(x) to convert between coordinates on the curve and the torus. The 
coefficients are meromorphic functions on 7~, and the corresponding Tracy-Widom kernel is 

K(XX) = nx)HY)-nYmx) (1124) 

A — Y 
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(ii) In [40] , Maier uses Lame's equation to produce explicit covering maps y — >■ T of the 
elliptic curve by hyperelliptic curves of arbitrary genus. For generic values of J and £ = 1,2,..., 
there exists a hyperelliptic curve 3^ of genus I and a holomorphic covering map 3^ — >■ T of 
degree + l)/2; in special cases, one can reduce hyperelliptic integrals to elliptic integrals. 

(iii) To solve the case I = 1, we introduce 

Mx, a) = -^((M^A), ^ - ^) JJ (1 _ g2 n )3) (n 25) 



n=l 



which satisfies Lame's equation — ^2^ 2 (^,a) + ^p(x)i/j2(x, a) = —p(a)ip2{x,a), and is such 
that a i ^ 4>2(x, a) is doubly periodic and x V^O^, a) is elliptic of the second kind; moreover 
■02 (#5 «)'02(— ol) = p(a) — p(x). By Lemma 11.2, fofa, oc) can be expressed as a quotient of 
tau functions from periodic linear systems as in Lemma 7.1, and Gaussian linear systems. 

(iii) (Integrable quantum systems) Having constructed the potential p from a periodic 
linear system, we can produce a family of Hankel kernels and potentials from standard lim- 
iting arguments which are associated with exactly solvable problems in quantum mechanics. 
Consider an interacting system of N identical particles at positions Xj on the real line which 
interact only pairwise, and where the strength of the mutual interaction of particles j and k 
depends only upon their separation Xj — Xk via a potential u; then the Hamiltonian is 

1 N d 2 

j=l 3 l<j<k<N 

Krichever [36] shows that the system with u(x) = 2p(x) is integrable in the sense that there 
exists a compact Riemann surface 3^v which covers the elliptic curve iV-fold, and the solution 
of the Hamiltonian dynamical system can be expressed in action-angle variables with the angles 
in the Jacobi variety of J^y. 

In each of the following, 7 and the potential u are meromorphic functions on a Riemann 
surface £ and ip satisfies the addition rule 

M x + y) = W x Mv\-M x WM . (H.27) 

7O) - i{y) 

£ u(x) 4>(x) j(x) t{x) 

P 1 g(g + l)/x 2 (g + l)/x -(g + l)/x 2 x a(9+i)/2 

C/nZ g(g + l)cosec 2 x (g + 1) cot x — (g + l)cosec 2 x (sinx) 9 ( 9+1 )/ 2 (11.28) 

C/niZ g(g + l)cosech 2 x (<7 + l)cothx — (g + l)cosech 2 x (smhx) g ^ 9+1 ^ 2 

C/A 2p(x \ A) ip 2 (x,a) -p(x \ A) d\(x | A) 

The rational, trigonometric and hyperbolic cases in this table satisfy the additional iden- 
tity u(x) = if^(x) 2 + 7(x) + c with c constant, and are derived from the quantum Lax equation 
in [p. 39, Su]. However, we can derive all of these as limiting cases from the elliptic potential 
in the final line of the table. We write A = Z2wi + Z2w 2 where u\,ui2/i > 0. Then we have 
the thermodynamic limit 

2p(x I Z2wi + Z2w 2 ) -» 2(7r/2w2) 2 cosech 2 (7r£/2cj 2 ) - 7r 2 /6u; 2 (wi -)• 00) (11.29) 
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and in contrast the high density limit 



2p(x | Z2a>i + Z2a>2) —> 2{tt/2ui\) cosec (ttx/2ui) - tt 2 /6ojI (uj 2 /i^oo); (11.30) 

in the latter case, q — >■ and hence 6\{x — a \ A)/8i(x | A) — > (sin7r(x — a))/ sin7rx, and 
one can check that the addition rule for ifj 2 reduces to the addition rule for the trigonometric 
function ip(x) = 2cot7rx. In each case, functions on T reduce to functions on the cylinder. 
When one limit is applied after the other, we have the limiting potential u(x) = 2/x 2 , and the 
addition rule gives the formula for Carleman's operator. See [48]. 

(iv) One can reverse the thermodynamic limit by taking a potential u to a periodic po- 
tential with period 2ui, as in u(x) t-> Xl^-oo u ( x + 2ooin). This observation is the basis for 
the scattering theory in [21, 32]. 

Any u G Kj- is associated with a finite collection of periodic linear systems as in Lemma 
11.2, and the spectrum of Schrodinger's equation is invariant under the KdV flow. We therefore 
consider how the motion of elliptic solitons under the KdV flow can be described in terms of 
periodic linear systems. 

Proposition 11.5 Let M be a finite dimensional differentiable manifold of elliptic functions 
on T that is invariant and differentiable with respect to the flow associated with KdV and 
that some u G M is finite gap. Then there exists a family E t = (—A, B(t), C; E(t)) of periodic 
linear systems such that u(x,t) is the potential from T, t , u(x,t) satishes KdV and S t evolves 
according to a hnite-dimensional Hamiltonian system. 

Proof. Airault, McKean and Moser showed [3, Corollary 1] that any such flow of potentials 
has the form 

m 

u(z,t) = ^2p(z- Xj{t)) +c (11.31) 
i=i 

where the moving poles Xj(t) lie on the manifold defined by the constraints 

m 

0= tp'ixj-xk) (fc = l,...,m). (11.32) 

and satisfy the system of nonlinear differential equations 

d m 

— =6 p(xj-Xk) (k = l,...,m) (11.33) 

In an evident analogy with (10.16), the Hamiltonian 

^ = ^5>f + y E p(*;-**) 2 ( n - 34 ) 

3 = 1 j^k:j,k = l 

gives this system of differential equations for the Xj\ see [17]. 
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We can realise 2p(x) as the potential of a periodic linear system (—A, B, C; E), and hence 
we can realise u(z, t) as the potential of the periodic linear system 



See [11] for more details of the construction and [17] for further information on the dynamics 
of the poles under KdV flows. 



Remark. We leave it as an open problem to characterize all finite gap cases of Hill's equation 
in terms of periodic linear systems. 
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